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Abstract 


A  first-order  detenalnatlon  of  the  effects  of  boundary  flow  of  high-speed 
aircraft  on  the  transmission  of  0.2^  to  20ij) radiation  is  presented.  The  four 
types  of  boundary  flow  considered  are:  main  flow,  shock  wave,  turbulent  boundary 
layer  (TBL),  and  separated  flow.  The  potential-flow  velocity  field  provides  a 
means  for  computing  wavefront  aberation  caused  by  curved  density  contours  in  the 
main-flow  region.  An  example  problem  la  worked  out  for  a  collimated  beam  emerging 
from  a  sphere.  The  density  jump  across  a  shock  wave  generally  causes  both  beam 
refraction  and  beam  degradation.  The  special  case  of  refraction  by  a  plane  shock 
wave  is  computed.  Estimates  of  the  average  modulation  transfer  function  (MTF)  for 
the  TBL  and  separated  flow  are  made  from  experimental  data.  From  the  wavefront 
aberation  of  main  flow  and  the  average  MTF' s  for  the  TBL  and  separated  flow,  the 
effect  of  boundary  flow  on  both  propagating  and  imaging  systems  Is  analyzed.  The 
effect  is  presented  in  terms  of  power  density  at  the  focus  point  for  propagating 
systems  and  in  terms  of  resolution  for  imaging  ay steps.  The  effect  of  shock  waves 
is  not  considered  in  the  analysis  because  the  only  quh^titativd  work  in  that  area 
concerns  beam  direction  only,  not  beam  quality.  lit  general  boundary  flow  has  only 
a  slight  effect  on  long  wavelength  (10.6)1)  radiation.  The  only  exception  is  the 
effect  of  main  flow  at  low  altitude  in  which  case  power  density  at  focus  may  drop 
to  less  than  one  percent  of  the  diffraction  limit.  A  similar  power  loss  may  be 
caused  by  both  main  flow  and  separated  flow  when  the  wavelength  is  short  (0.55)1). 
Only  the  TBL  effect  is  not  drastic.,.  Finally,  the  scaling  rules  for  predicting 
boundary-flow  optical  effects  of  a  Body  from  a  similar  body  of  different  scale  is 
derived.  Adjustment  of  wavelength  and/or  density  is  required  to  compensate  for  the 
change  in  scale. 
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I.  Introduction 


The  quality  of  an  optical  beam  which  passes  through  the  boundary  flow  of  a 
transonic  or  supersonic  aircraft  can  be  drastically  reduced.  Camera  resolution 
dropped  from  164  lines/mm  to  30  lines/mm  in  a  wind  tunnel  test  of  the  RF-4C 
aircraft.1  When  a  lower  acuity  system  was  flight  tested  in  the  RF-4C,  resolution 
dropped  from  35  lines/mm  to  20  lines/mm. 2  The  reason  for  the  resolution  loss  was 
that  the  optical  path  from  the  target  passed  through  a  separated  flow  region 
behind  the  nose  chin  as  shown  in  Figure  1.  Density  fluctuations  in  this  region 
caused  random  wavefront  errors,  and  resolution  was  thus  reduced.  When  the 
aircraft  nose  was  refaired,  the  resolution  loses  did  not  occur. 3,4  Random 
density  fluctuations  are  not  the  only  source  of  aerodynamically  induced  optical 
degradation.  It  will  be  shown  that  the  curved  air  density  contours  in  the  main 
flow  region  of  a  blunt  body,  as  shown  in  Figure  2,  produce  wavefront  aberrations 
which,  if  uncorrected,  may  completely  destroy  the  optical  quality  of  an  incoming 
or  propagated  beam. 


These  examples  show  that  an  understanding  of  boundary  flow  effects  is 
necessary  if  airborne  optical  systems  are  to  achieve  optimal  results. 

Purpose  of  this  Report 

The  purpose  of  this  report  is  to  present  a  first-order  determination  of  the 
effects  of  boundary  flow  on  the  transmission  of  0.2u  to  20p  radiation.  This 
abroach  establishes  the  significance  of  the  aerodynamic  effect  and  hopefully  will 
lead  to  more  rigorous  analytical  and  experimental  studies. 

Types  of  Boundary  Flow 

Boundary  flow  may  be  divided  into  two  classes,  random  and  controlled.  In 
addition  to  separated  flow,  random  flew  includes  the  turbulent  boundary  layer  (TBL) , 
which  forms  on  the  skin  of  all  aircraft.  The  TBL  causes  a  degradation  similar  to, 
but  weaker  than  that  caused  by  separated  flow.  Stine  and  Vinovich  have  shown  that 
TBL's  as  thin  as  one  inch  may  cause  significant  degradation.5  Unlike  the  random 
flow,  controlled  flow  has  density  gradients  which  are  constant  in  time.  Also 
wavefront  deviations  caused  by  controlled  flow  generally  have  a  larger  spatial 
scale  than  those  caused  by  random  flow  as  shown  in  Figure  3.  Both  main  flow  and 
shock  waves  are  types  of  controlled  flow.  The  shock  wave  occurs  whenever  some  part 
of  the  flow  field  outside  an  aircraft  has  a  speed  greater  than  the  speed  of  sound 
relative  to  the  aircraft.  The  shock  wave  is  a  very  thin  pressure  jump  and  there¬ 
fore  a  very  thin  density  jump.  In  the  simplest  kind  of  wavefront  deviation,  the 
shock  wave  is  a  plane  interface  between  media  of  different  refractive  indices,  and 
a  beam  passing  through  the  shock  will  behave  according  to  Snell's  law  as  shown  in 
Figure  4.  All  four  types  of  boundary  flow  (separated  flow,  TBL,  main  flow,  and 
shock  wave)  are  shown  in  Figure  5  for  an  RF-4C  at  Mach  .93. 
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FIGURE  1 

WAVEFRONT  DEGRADATION  CAUSED  BY  FLOW  SEPARATION 
BEHIND  RF-4C  NOSE  CHIN 
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FIGURE  2 

WAVEFRONT  DEGRADATION  CAUSED  BY  CURVED  AIR-DENSITY  CONTOURS  IN 
THE  MAIN-FLOW  REGION  OF  A  BLUNT  BODY 
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FIGURE  3 

COMPARISON  OF  CONTROLLED- AND  RANDOM-WAVE  FRONT  DEGRADATION 
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FIGURE  5 

POSITIONS  OF  MAIN  FLOW,  SHOCK  WAVE,  TURBULENT  BOUNDARY  LAYER,  AND 
SEPARATED  FLOW  FOR  AN  RF-4C  AT  MACH  0.93 
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This  report  considers  the  optical  effects  of  the  four  types  of  boundary  flow: 
main  flow,  shock  wave,  TBL,  and  separated  flow.  The  next  section  presents  the 
basic  aerodynamics  and  optics  theory  used  in  the  rest  of  the  report. 

Controlled  wavefront  deviations  are  discussed  in  Section  III.  Wavefront 
aberration  cuased  by  curved  density  contours  in  the  main-flow  region  are  computed 
from  the  potential-flow  velocity  field.  An  example  problem  is  worked  out  for  a 
collimated  beam  emerging  from  a  sphere.  This  section  also  includes  a  qualitative 
description  of  degradation  caused  by  shock  waves.  The  special  case  of  refraction 
by  a  plane  shock  wave  is  computed. 

In  Section  IV  the  Stine-Winovich  experiment  data  and  the  RF-4C  wind  tunnel 
data  are  used  to  estimate  average  modulation  transfer  functions  (MTF's)  for  TBL's 
and  separated  flow  respectively.  A  recent  experiment  by  Kelsall  is  also  described. 
Althoug’  "iveral  complications  make  it  impossible  to  properly  analyze  the  resultant 
data,  this  experiment  is  important  because  the  average  boundary  flow  MTF  was 
measured  directly  using  a  shearing  interferometer.  Hopefully  this  new  technique, 
which  produces  more  reliable  MTF  measurements  than  any  other  method,  will  be  used 
in  controlled  experiments  to  measure  boundary  flow  MTF's. 

From  the  results  of  the  previous  sections,  Section  V  analyzes  the  optical 
effect  of  boundary  flow  on  both  propagating  and  imaging  systems.  The  effect  is 
presented  in  terms  of  power  density  at  the  focus  point  for  propagating  systems  and 
in  terms  of  resolution  for  imaging  systems.  In  addition  to  main  flow,  turbulent 
boundary  layers,  and  separated  flow,  wake  flow  is  also  considered.  Wake  flow  is 
simply  separated  flow  which  occurs  behind  a  blunt  body  as  shown  in  Figure  6.  In 
Section  V  wake  flow  is  treated  as  a  special  case  of  separated  flow.  Shock  waves 
are  not  considered  because  the  only  quantitative  work  done  in  that  area  concerns 
beam  direction  only,  not  beam  quality.  Finally,  Section  V  derives  the  scaling  rules 
for  predicting  boundary-flow  optical  effects  of  a  body  from  a  similar  body  of 
different  scale. 

Notation,  Symbols,  and  Important  Equations 

Appendix  A  contains  a  summary  of  important  equations  as  well  as  a  list  of 
notations  and  symbols. 
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II.  Basic  Theory 

The  purpose  of  this  section  is  to  present  the  basic  aerodynamic  and  optics 
theory  used  in  the  rest  of  the  report.  The  topics  discussed  here  are  the 
relationship  between  refractive  index  and  density  for  air,  the  optical  phase,  the 
optical  transfer  function  (OTF),  power  density  in  the  focal  plane,  and  the  aero¬ 
dynamic  similarity  rules. 

Relationship  Between  Refractive  Index  and  Density  for  Air 

Because  the  refractive  index  n  of  air  is  very  close  to  unity,  it  is  conven¬ 
ient  to  define  the  refractive  modulus 

N  -  (n-1)  x  106.  (2-1) 

At  normal  temperatures  and  pressures  N  is  of  the  order  of  300.  A  formula  which 
gives  N  to  within  +  1  for  dry  air  in  the  0.2u  to  20p  wavelength  region  is 

N  -  0.776  f  (1  +  7,^2-10—5)  (2-2) 

™  A 


where 

2 

P  ■  pressure  in  nt/m 
T  -  temperature  in  °K  7 
X  -  wavelength  in  meters. 

A  correction  factor  for  water  vapor,  which  may  subtract  as  much  as  3  from  N,  is 
given  by 

-  0.6  x  10'3  Pw 

where  P  is  the  water  vapor  pressure  in  nt/m^.  For  the  accuracy  required  in  this 
paper,  Kowever ,  the  correction  factor  will  not  be  required.  Using  the  perfect  gas 
law, 

P  -  287  pT  (2-3) 

where  p  is  density  in  kg/m^,  Equation  (2-2)  may  be  written 

N  -  223  P  (1  +  -,52  *  —  15 )  (2-4) 

X* 

The  refractive  index  may  now  be  written 

n  -  1  +  Kp  (2-5) 


where  the  Glads tone-Dale  constant  K  is  given  by 


K  -  0.223  x  10"3  (1  + 


7.52  x  10 
X2 


-15 


)  m3/kg 


(2-6) 
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The  wavelength  term  is  negligible  for  wavelengths  longer  than  10W  so  that 

K  -  0.223  x  l(f3  m3/kg  (2-7) 

for  long  wavelengths.  Even  for  shorter  wavelengths,  the  wavelength  term  may 
have  little  effect;  for  example,  when  X  ■  0.55y  K  increases  to  only  0.228  x  10 
nrtykg.  For  most  purposes,  therefore,  Equations  (2-5)  and  (2-7)  may  be  used  to 
calculate  n.  Otherwise  Equations  (2-5)  and  (2-6)  may  be  used. 

Optical  Phase 

For  our  purposes  we  may  represent  a  general  optical  wave  as  a  scalar  given 
by 

F(x,y,z,t)  ■  A(x,y,z,t)  exp{-i[<ut  -  kA(x,y,z,t)] }  (2-8) 

where 

x,y,z  ■  Cartesian  coordinates  in  space 
t  ■  time 

F  *  wave  function  (generally  complex) 

A  -  wave  amplitude  (generally  complex) 

A  ■  wave  phase 

u  "  angular  frequency  (2rrc/X) 
k  -  wave  number  (2ir/X) 
c  “  speed  of  light 
X  -  optical  wavelength. 

The  function  F  gives  a  complete  description  of  the  optical  wave  in  space  and  time. 
In  this  subsection  we  will  show  that  for  our  case  the  quotient  A/X  completely 
describes  the  optical  wave.  The  relationship  between  A  and  the  refractive  index  n 
of  the  medium  is  given.  Finally,  it  will  be  shown  that  A  and  the  wavefront  W  are 
equal  in  magnitude  and  opposite  in  sign. 

First  Equation  (2-8)  is  simplified.  Since  our  medium  is  not  lossy  and  since 
the  wavefront  is  essentially  plane,  not  spherical,  A  is  a  constant.  The  term 

-iwt 

e 

in  Equation  (2-8)  is  left  to  be  understood.  The  time  dependence  of  A  is  also  left 
to  be  understood,  and  Equation  (2-8)  becomes 

FC.y.z)  -  (2-9) 

This  equation  shows  that  at  a  fixed  time  and  amplitude  F  depends  on  A/X  only. 

Let  us  examine  a  beam  which  traverses  a  medium  of  varying  refractive  index. 

A  plane  parallel,  monochromatic  beam  emerges  from  the  x'y'-plane  of  a  Cartesian 
coordinate  system  with  its  direction  of  propagation  parallel  to  the  z'-axis  as 
shown  in  Figure  7.  At  some  distance  L>>X  from  the  x'y '-plane  we  examine  the 
wave. 


REPORT  MDC  A2582 
14  DECEMBER  1973 


x' 


BEAM  TRAVERSING  MEDIUM  OF  VARYING  REFRACTIVE  INDEX 


REPORT  MOC  A2582 
14  DECEMBER  1973 


F(x,y,L)  -  AelkA<X,y,L>  (2-10) 

The  phase  is  found  by  integrating  the  refractive  index  n  over  the  optical 

path. 

The  phase  is  given  by 


(x.£,L) 

A(x,y,L)  *  j  n(x',y',z')  ds  +  a  X  (2-11) 

( x,$,0 ) 

where  ds  is  an  increment  along  the  ray  path  from  (xty,0)  to  (x,y,L)  and  m  is  any 
convenient  integer.  For  the  boundary  flow  problem  the  refractive  index  changes 
are  small,  so  the  ray  path  will  be  very  close  to  parallel  to  the  z'-axis.  The 
variations  of  the  ray  from  the  z'-axis  in  Figure  7  are  greatly  exagerated. 

ds  -  dz'  (2-12) 


and 


A(x,y,L) 


n(x,y,z')  dz'  +  m  X. 


(2-13) 


We  now  look  at  the  wavefront  W(x,y,L)  at  z'  ■  L  as  shown  in  Figure  7.  It  is 
defined  by  specifying  that  the  phase  at  the  wavefront  is  equal  to  some  arbitrary 
constant. 


W(x,y,L) 

A[x,y,W(x,y,L) ]  -  const  -  f  n(x,y,z')  dz'  +  m  X  (2-14) 

or 

L  W(x.y,L) 

const  ■  I  n(x,y,z')  dz'  +  m  X  +  f  n(x,y,z')  dz'  (2-15) 

Using  Equation  (2-13)  and  the  fact  that  n  is  constant  over  distances  of  several 
wavelengths.  Equation  (2-15)  becomes 

const  -  A(x,y,L)  +  n(x,y,L)  W(x,y,L)  (2-16) 

Choosing  zero  for  the  value  of  the  arbitrary  constant  and  noting  that  n  is  very 
close  to  unity,  we  have 

W(x,y,L)  -  -A(x,y,L) .  (2-17) 

Equation  (2-17)  shows  that  4  and  W  have  the  same  magnitude  but  opposite  signs. 
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Optical  Transfer  Function 


In  the  case  of  turbulent  boundary  layers  and  separated  flow,  A  is  a  random 
variable  which  ifl  practice  cannot  be  measured.  However,  the  average  OTF,  which 
is  a  function  of  A,  can  often  be  measured  or  computed. 


Definition.  If  in  Figure  7  a  perfect  lens  is  placed  perpendicular  to  the 
z'-axis  at  z'  -  L,  an  image  is  formed  in  the  focal  plane.  A  screen  with  an 
aperture  £  is  placed  adjacent  to  the  lens  as  shown  in  Figure  8  so  that  all  light 
passing  through  the  aperture  must  also  pass  through  the  lens.  The  pupil  function 
G  is  defined  in  terms  of  the  wave  F  at  the  aperture  as  follows. 


G(x,y) 


F(x,y)  -  A.iki<*’y),  (x,y)e£ 
|  o,  (x,y)^  r 


(2-18) 


In  words,  the  pupil  function  G  is  identical  to  the  wave  function  F  inside  the 
aperture  but  vanishes  outside  the  aperture.  The  dependence  of  F  and  A  on  L  is 
left  to  be  understood.  The  wave  amplitude  U  in  the  focal  plane  is  the  inverse  Fourier 
transform  of  G, 


U(X,Y)  -  ^  J j"  G(x,y)  exp£-^(xX  +  yY)]  dx  dy  (2_19) 

where  X  and  Y  are  the  coordinates  in  the  focal  plane  and  f  is  the  focal  length. 

The  image  irradiance  X  is  then  given  by 

I (X,Y)  -  U*(X,Y)  U(X,Y) .  (2-20) 

The  asterisk  denotes  complex  conjugate.  The  optical  transfer  function  t  of  the 
wave  passing  through  the  aperture  is  defined  to  be  the  normalized  Fourier  transform 
of  I, 


'(»•»)  ■  ?  fj  I(X'Y>  + 


dX  dY 


(2-21) 


where 


00 

■Si 


I(X,Y)  dX  dY 


(2-22) 


Using  Parseval's  theorem  from  Fourier  analysis  in  conjunction  with  Equations  (2-19), 
(2-20),  and  (2-22)  we  have 


ii 


U*(X,Y)  U(X,Y)  dX  dY 


■A 


G*(x,y)  G(x,y)  dx  dy 


(2-23) 
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FIGURE  8 

GEOMETRY  FOR  DEFINING  OPTICAL  TRANSFER  FUNCTION 
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or  using  Equation  (2-18) ,  we  have 


C  -  A*AE . 


(2-24) 


Equations  (2-22)  through  (2-24)  each  imply  that  C  is  the  radiant  power  which  passes 
through the  aperture  and  lens  to  the  £ocal  plane. 

The  spatial  frequencies  defined  by 


vx  "  If  and  vy  “  Xf 


(2-25) 


allow  Equation  (2-21)  to  be  written 


r(Vx'Vy)“  C  JJI(X*Y)  exp[±2,r(] 


Xv  +  Yv  I  dX  dY. 

l  x  y/J 


(2-26) 


The  quantities  v  and  vy  have  units  of  cycles  per  unit  distance.  Since  t  is  the 
Fourier  transform  of  a  real  function,  it  is  Kermition;  the  real  part  of  t  is  an 
even  function,  and  the  imaginary  part  is  odd.  If  I  is  even,  r  is  real.  The 
MTF  can  now  be  defined  simply  as  the  absolute  value  of  the  OTF. 

Autocorrelation  Function  of  Pupil  Function.  Although  the  OTF  is  defined  by 
Equation  (2-21),  it  is  often  convenient  to  express  x  as  a  function  of  the  pupil 
function  rather  than  the  focal  plane  intensity.  Specifically,  it  can  be  shown 
that  the  OTF  is  equal  to  the  normalized  autocorrelation  function  of  the  pupil 
function  by  using  a  theorem  from  Fourier  analysis. 


^TPX?)  "  C  JJ  W*'**')  G(x'+x,y'+y)  dx'dy.' 


(2-27) 


The  OTF  may  be  expressed  in  terms  of  A  by  using  Equations  (2-18)  and  (2-24), 


■(*$)  ■  r  JJ 


GQ(x',y')  GQ(x'+x,y'+y) 


x  exp{ik[A(x'+x,y'+y)  -  A(x',y')]}  dx'dy' 


(2-28) 


where 


GQ(x,y) 


1.  (x,y)  t  E 


0,  (x.y)  i  I 


(2-29) 


Diffraction  OTF.  A  special  case  occurs  when  the  wave  front  entering  the  aper¬ 
ture  is  plane,  in  which  case  A  is  constant.  Equation  (2-28)  becomes  the  diffraction- 
limit  OTF  (or  MTF,  since  the  OTF  is  real)  of  the  aperture. 
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To(t?*Xf)  "  2  G0(x'*y>)  G0<*'+*'j'+y)  to'dy'  (2-30) 

Two  aperture  shapes  which  are  often  used  are  the  rectangle  and  the  circle.  For 
a  retangular  aperture  with  dimensions  Dx  x  the  pupil  function  is 


GQ(x,y) 


1»0<_x<.Dx,  0  <.  y  <.  Dy 
0  otherwise 


(2-31) 


and  the  diffraction  OTF  is 


-D  <  x  <  D  ,  -D  <y<D 
x—  —  x  y  —  —  y 

(2-32) 


For  a  circular  aperture  of  diameter  D  the  pupil  function  is 


G0(D 


1,  0  <_  r  <_  D/2 
0  otherwise 


(2-34) 


where 


r  ■  *^x^  +  y^. 


(2-34) 


The  diffraction  OTF  is 


T 


0 


0  <  r  <  D 


(2-35) 


Figure  9  is  a  plot  of  the  diffraction  OTF's  for  rectangular  and  circular  apertures. 
Notice  that  only  positive  spatial  frequencies  are  used  since  all  diffraction  OTF’s 
are  even  functions. 


Examination  of  Equation  (2-28)  reveals  that  any  aberation  in  the  wavefront 
at  the  aperture  can  only  serve  to  decrease  the  MTF  below  the  diffraction  value. 
That  is 


(2-36) 


A  comparison  of  t  with  t  is  often  used  as  a  measure  of  the  seriousness  of  the 
aberation. 


Diffraction  OTF  7 
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FIGURE  9 

DIFFRACTION  OTF  FOR  RECTANGULAR  AND  CIRCULAR  APERTURES 

GP7*-0?0O-9 
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Sine-Wave  Response  Function.  For  an  imaging  system  the  MTF  is  the  sine-wave 
response  function  of  the  system.  Consider  the  Imaging  system  geometry  shown  in 
Figure  10.  An  incremental  area  dX'dY"  located  at  (X'.Y')  in  the  object  plane 
produces  an  image  on  the  image  plane.  We  assume  that  light  from  the  object  plane 
is  incoherent,  and  for  simplicity  we  assume  that  the  magnification  is  unity.  The 
radiant  power  reaching  the  image  plane  from  our  incremental  area  is  proportional 
to  0(X',Y')  dX'dY'  where  0(X',Y')  is  the  object  radiance  at  (X",Y').  We  now  wish 
to  find  the  resulting  image  irradiance  at  (X,Y).  In  defining  the  OTF,  I(X,Y)  was 
the  image  irradiance  resulting  from  a  wave  which,  except  for  boundary  flow  effects, 
was  plane  at  the  lens  aperture.  We  may  also  take  I(X,Y)  as  the  image  irradiance 
resulting  from  a  point  object  far  to  the  left  of  the  lens.  In  this  sense  we 
now  define  the  point  spread  function  p(X,Y)  as  the  normalized  image  irradiance 
corresponding  to  a  point  object.  The  OTF  is  now  taken  as  the  Fourier  transform 
of  p(X,Y).  In  the  present  context  I(X,Y)  results,  not  from  a  point,  but  from  a 
finite  object.  It  is  this  image  irradiance  which  we  are  trying  to  find.  Since 
the  incremental  area  may  be  considered  a  point,  the  incremental  contribution  to 
I (X,Y)  is 

dl (X,Y)  -  0(X",Y')  dXMY'  p(X-X'.Y-Y') .  (2-37) 

We  assume  that  p(X,Y)  is  constant  for  object  points  throughout  the  object  plane. 

I (X,Y)  is  found  by  integrating  Equation  (2-37)  over  the  object  plane. 


00 

I (X,Y)  -  JJ  O(X'.Y')  p(X-X'.Y-Y')  dX'dY' 


(2-38) 


The  integral  is  called  the  convolution  of  0  and  p.  By  applying  the  convolution 
theorem  from  Fourier  analysis  to  Equation  (2-38) Jwe  have 


3.  (v  ,v  J  *  S  (v  ,v  }  x  L  ,v  } 

i\x  y)  °\*  yj  \x  yJ 


(2-39) 


where  S^,  S  ,  and  t  are  the  Fourier  transforms  of  I,  0,  and  p  respectively. 
Equation  (2-39)  shows  that  the  OTF  is  the  spatial  frequency  response  of  the 
system. 

We  now  let  the  object  be  a  sinusoid  of  frequency  v, 

0(X',Y')  *  C0  sin(2wvX'+  *)  +  1,  0  <  C  <  1 


(2-40 


where  4>  is  a  phase  factor  and  Cq  is  a  constant  which  we  call  the  object  modulation. 
Taking  the  Fourier  transform  we  have 


s»(v°)  ■  s  *{\  - ») + i\) 


(2-41) 
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IMAGING  SYSTEM  GEOMETRY 
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where  <$  Is  Che  delta  function.  Substituting  Equation  (2-41)  into  Equation  (2-39) , 
we  have 


Si(V°)!  ■  I  t(vx,0)  I c0  5(vx-v)  +|t(vx,0)|5(vx) 
■  | t(v,0) |Cq6(vx-v)  +  fi(vx). 


(2-42) 


This  equation  shows  that  the  image  of  a  sinusoid  of  frequency  v  is  also  a  sinusoid 
of  frequency  v,  the  modulation  of  which  is  reduced  by  a  factor  of  |t(v,0)|.  In 
other  words,  the  MTF  is  the  sine-wave  response  function  of  an  imaging  system  as 
illustrated  In  Figure  11. 

Power  Density  in  Focal  Plane 

Suppose  a  beam  of  light  emerging  from  and  aircraft  is  focussed  at  some  distant 
point  as  shown  in  Figure  12.  Diffraction  and  the  boundary  flow  will  cause  the 
point  image  to  be  spread  over  the  focal  plane.  In  this  type  of  problem  the  image 
irradiance  or  power  density  I(X,Y)  is  sought.  The  problem  geometry  here  is 
somewhat  different  than  that  used  earlier.  Here  the  optical  distrubance  caused 
by  boundary  flow  occurs  between  the  lens  and  the  focus  point. 

The  mathematics  developed  for  the  geometry  in  which  the  lens  is  between  the 
disturbance  and  the  focus  point  is  still  valid,  however,  as  long  as  the  distance 
between  the  lens  and  the  disturbance  is  much  less  than  the  focal  length. 

One  way  to  attack  this  problem  would  be  to  use  Equations  (2-19)  and  (2-20) . 
Unfortunately  the  pupil  function  G  is  often  unknown.  Since  the  average  OTF  may 
be  known,  however,  we  look  for  a  way  to  compute  I  from  t.  This  is  easily  found 
since  Equation  (2-21)  shows  that  I  and  x  are  Fourier  transform  pairs. 


I(X,Y) 


(Xf)‘ 


[{  T(n’&)exp['  f (XX  +  yY)] 


dx  dy  (2-43) 


In  polar  coordinates  with  rotational  symmetry  we  have 

«*>  '  (ITT?  |  '(*)  >o{¥) ' 

where 


R 


/x2  +  y2 


(2-45) 
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D  is  the  aperture  diameter,  and  J.  Is  the  zero  order  Bessel  function  of  the  first 
kind.  Since  Integration  is  a  linear  process,  t  and  I  may  be  replaced  by  their 
averages  in  Equations  (2-43)  and  (2-44). 

Often  the  average  power  density  at  the  beam  center  is  of  particular  interest. 
From  Equations  (2-43)  and  (2-44)  we  have 


1(0,0)  - 


(Xf)‘ 


5 


dx  dy 


and 


1(0)  - 


2ffC 

(Xf)5 


H*) 


r  dr . 


(2-46) 


(2-47) 


In  Section  V  we  use  the  ratio  of  1(0,0)  to  the  central  power  density  achieved 
in  the  diffraction  limit,  I-(0,0).  By  substituting  Equation  (2-30)  into  Equation 
(2-46) ,  we  find  that 


I„(0,0) 


Cl 


(Xf) 


2  • 


The  ratio  is  then 


KO.O)  Iff  ^ 

yo!o)  i  )3  ^xrxf  jdx  ** 


(2-48) 


(2-48) 


or 


K0) 

v°) 


D2  I  '(») 


r  dr. 


(2-50) 


In  this  report  the  ratio  above  is  called  the  normalized  central  power  density 
(NCPD) . 

The  NCPD  can  also  be  expressed  in  terms  of  the  optical  phase  A  by  using 
Equations  (2-18),  (2-19),  (2-20),  (2-24)  and  (2-48). 


NCPD 


Sj  V*.y> 


(2-51) 


If  we  have  polar  coordinates  with  rotational  symmetry,  Equation  (2-51)  becomes 

D/2 


NCPD  -  J  eikA(r)r  dr|2. 


(2-52) 


1 


MST 
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Aerodynamic  Similarity  Rule 

£ 

The  aerodynamic  similarity  rule,  which  is  taken  from  fluid  flow  theory, 
states  that  two  similar  models  of  different  scale  will  produce  similar  flow  if 
certain  conditions  are  met. 9  This  rule  is  used  in  developing  optical  scaling  rules 
in  Section  V. 

In  general,  the  fluids  for  the  two  models  are  different,  and  there  are  six 
conditions  that  must  be  met.  Before  these  conditions  are  presented  it  is  necessary 
to  introduce  some  notation.  Let  us  define  the  following  terms  for  any  model. 

X  characteristic  linear  dimension  (for  instance,  the  distance  between  two 
given  points  on  the  model) 

t  characteristic  time  (for  example,  the  time  it  takes  for  a  fluid  molecule 
to  travel  along  the  model  surface  from  one  given  point  to  another) 

V  fluid  speed 

cp  fluid  specific  heat  at  constant  pressure 

cv  fluid  specific  heat  at  constant  volume 

a  speed  of  sound  in  the  fluid 

g  gravitational  acceleration 
p  fluid  density 

p  fluid  coefficient  of  viscoisty 

<  fluid  coefficient  of  thermal  conductivity 

General  Conditions .  Two  similar  models  of  different  scale  will  produce  similar 
flows  if  the  following  conditions  are  met. 

1.  Corresponding  values  of  the  time  scale  (X/tV)  are  equal 

2.  Corresponding  values  of  the  ratio  of  specific  heats  (y-Cp/ cv)  are  equal 

3.  Corresponding  values  of  the  mach  number  (M“V/a)  are  equal 

4.  When  gravitational  body  forces  on  the  fluid  are  important,  corresponding 
values  of  the  Froede  number  (gX/V?)  are  equal. 

5.  Corresponding  values  of  the  Reynolds  number  (R*  ■  pVX/u)  are  equal 

6.  Corresponding  values  of  the  Prandtl  number  (ucp/ic)  are  equal. 

Although  there  are  six  conditions  which  generally  must  be  met,  only  three  are 
important  with  respect  to  high  speed  aircraft. 

Pertinent  Conditions.  Because  the  change  in  gravitational  acceleration  is 
negligible  over  the  entire  region  near  the  aircraft,  gravitational  body  forces  on 
the  fluid  are  not  important  and  Condition  4  can  be  neglected.  Condition  2  can  be 
neglected  because  the  fluid  under  consideration  is  always  air  which  has  a  constant 


*  We  refer  here  to  the  differential  equations  of  motion  of  a  viscous,  compressible, 
heat  conducting  fluid.  Because  some  of  these  equations  are  nonlinear  partial  dif¬ 
ferential  equations,  a  general  solution  has  never  been  found.  The  interested 
reader  is  referred  to  a  text  on  fluid  mechanics  such  as  References  9  and  10. 
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Y  of  1.4.  Finally  Condition  6  may  be  neglected  because  a  constant  Y  Implies  a 
constant  Prandtl  number.’  The  remaining  Conditions  1,  3,  and  5  which  are  pertinent 
to  high  speed  aircraft  state  that  corresponding  values  of 

Time  scale  (X/tV) 

Mach  Number  (V/a) 

Reynolds  number  (pVX/y) 

are  equal  for  the  two  models.  The  values  of  a  and  y,  which  are  functions  of 
temperature  only,  are  given  by 


/287  yT  m/sec 


(2-53) 


and 


where  T  is  temperature  in  °K. 
1,  3,  and  5  we  require 


kg/m-sec 


(2-54) 


Consider  now  two  models,  1  and  2.  To  meet  Conditions 


X1/C1V1  *  X2ft2V2 

(2-55) 

Vl/al  *  V2/a2 

(2-56) 

plVlXl/ul  "  P2V2X2/U2 

(2-57) 

If  the  temperatures  for  the  two  models  are  equal,  then 
Equations  (2-55),  (2-56),  and  (2-57)  imply 

ai  -  a2,  Pj  ■  y2,  and 

V2  -  Vx  or  Mj_  -  M2 

(2-58) 

X2 

t2  “  C1  X1 

(2-59) 

and 

X1 

P2  *  P1  X2 

(2-60) 

This  means  that  if  we  wish  a  one-fifth  scale  model  (Model  2)  to  produce  a  flow 
similar  to  that  of  a  full  size  airplane  (Model  1) ,  then  the  speed  must  be  kept  the 
same,  and  the  density  must  be  increased  by  a  factor  of  five.  The  characteristic 
time  of  the  flow  will  decrease  by  a  factor  of  five. 
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III.  Controlled  Effects 

This  section  presents  a  discussion  of  the  controlled  wavefront  deviations 
caused  by  the  main-flow  region  and  by  shock  waves. 

Main  Flow 

This  subsection  presents  a  general  method  for  predicting  the  optical  phase 
disturbance  A(x,y)  caused  by  density  variation  in  the  main-flow  region  of  a  sub¬ 
sonic  body.  This  method  is  then  applied  to  the  main -flow  field  about  a  sphere. 

General  Method.  Figure  13(a)  demonstrates  the  problem.  A  beam  emerges 
from  an  aircraft,  traverses  the  main-flow  region,  and  is  focused  at  a  point 
very  far  away.  The  wavefront  emerging  from  the  aperture  is  assumed  to  be  plane 
parallel.  He  wish  to  predict  the  wave  phase  A(x,y)  at  a  distance  L  from  the 
aircraft  such  that  density  beyond  L  is  essentially  equal  to  freestream  density. 
The  solution  to  this  problem  is  also  the  solution  to  the  problem  represented  by 
Figure  13(b).  We  are  given  a  plane  parallel  wavefront  at  L  moving  toward  the 
aircraft,  and  we  wish  to  know  A(x,y)  at  the  aircraft. 

We  begin  by  stipulating  that 


A(0,0)  -  0  (3-1) 

where  the  origin  is  taken  as  the  center  of  the  outgoing  wavefront.  The  phase  may 
then  be  expressed  as 


L 

M*.y)  -  $  [n(x,y ,z)  -  n(0,0,z)]  dz  (3-2) 

0 

where  n  is  refractive  index.  Substituting  Equation  (2-5)  into  Equation  (3-2),  we 
have 


L 

A(x,y)  -  K  J  [p(x,y,z)  -  p(0,0,z)  ]  dz  (3-3) 

0 

where  p  is  density  and  K  is  a  constant  given  by  Equation  (2-6)  or  approximately  by 
Equation  (2-7).  In  order  to  solve  Equation  (3-3)  it  is  necessary  to  find  the 
density  field  about  the  aircraft  (neglecting  the  boundary-layer  and  separated-flow 
regions).  This  can  be  done  by  first  finding  the  velocity  field  for  an  incompres¬ 
sible  fluid  using  potential  flow  theory  and  then  applying  the  Karman-Tsien  hodograph 
method  as  presented  by  Liepman  and  Puckett. H 
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In  working  this  type  of  problem  it  is  helpful  to  use  the  reference  system 
shown  in  Figure  14.  The  aircraft  is  stationary  as  is  a  very  large  fictitious 
reservoir  of  compressed  air.  A  large  hole  in  the  reservoir  allows  air  to  escape 
and  pass  by  the  aircraft  at  a  free  stream  speed  U.  We  let  ag  and  aa  represent  the 
speed  of  sound  in  the  reservoir  and  in  the  frees tream  respectively,  and  we  let 
Pq  and  represent  the  air  density  in  the  reservoir  and  in  the  freestream 
respectively.  The  quantities  U,  aB,  and  p^  are  given. 

Using  the  shape  of  the  aircraft  and  the  given  value  of  U,  the  velocity  square 
field  v2(x,y,z)  for  an  incompressible  fluid  is  found  from  potential  flow  theory. 
General  methods  for  calculating  v2  are  not  considered  here,  except  for  the  following 
comment.  The  field  can  be  easily  found  for  some  simple  shapes  such  as  for  a 
sphere.  For  more  complicated  shapes,  solutions  can  be  found  by  superimposing  the 
velocity  potentials  of  simple  shapes. ^ 

Once  V2(x,y,z)  is  found,  the  density  field  p(x,y,z)  can  be  found  by  the 
following  approximation  from  the  Karman-Tslen  hodograph  method. 11 

4a3  -  V2 (x,y,z) 

p(x,y,z)  ■  P0  — 5 - z -  (3-4) 

4ag  +  VZ(x,y,z) 

This  equation  may  be  rewritten  in  terms  of  pw  and  a*  instead  of  pQ  and  ag  by 
applying  the  equalities11 


where  y  is  the  ratio  of  specific  heats.  For  air  y  “  1.4.  Using  Equations  (3-4), 
0-5),  and  (3-6),  Equation  (3-3)  may  be  expressed  as 


M*,y>  -  h  J 


w'(x,y,z) 


1  -  w'(0,0,z)  iz 
1  +  w'(0,0,z) 


(3-7) 


where 

H  -  Kp,  U  +  ^  Mfl. 


(3-8) 


M  “  U/a  “  freestream  Mach  number. 


(3-9) 
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Speed  of  Air  =  0 
Speed  of  Sound  *  ae 
Air  Density  ■  pQ 


Aircraft 


Free-Straam  Speed  of  Air  *  U 
Spaed  of  Sound  «  a„ 

Air  Density  «  pQ 


FIGURE  14 

REFERENCE  SYSTEM  FOR  MAIN-FLOW  PROBLEM 
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and 


w'(x,y,z)  - 


V2(x,y,2)/af 
2(y-1)M^  +  4 


(3-10) 


Equations  (3-7)  through  (3-10)  provide  a  method  for  estimating  A(x,y)  provided 

V2(x,y,z)  <  a2  (3-11) 

for  all  space.  The  speed  aB  is  given  by10 


(3-12) 


where  T^  is  the  frees  treason  temperature  in  8K  and  a,,,  is  in  meters /sec. 

Main  Flow  for  Sphere.  To  find  A(x,y)  for  a  sphere  the  coordinate  system  shown 
in  Figure  15  is  used.  The  polar  axis  points  into  the  flow,  and  the  center  of  the 
beam  defines  the  z-axis.  The  angle  between  the  polar  axis  and  the  z-axis  is  0g. 

The  y-axis  lies  perpendicular  to  the  z-axis  in  the  plane  defined  by  the  polar  and 
z  axes.  The  x-axis  is  perpendicular  to  both  the  y  and  z  axes.  The  beam  emerging 
from  the  sphere  is  taken  to  be  plane  parallel,  lying  paralled  to  the  x,y-plane, 
at  z-R  where  R  is  the  sphere  radius.  To  simplify  the  problem  the  following 
normalized  variables  are  Introduced. 


C  -  x/R  t\  -  y/R  S  -  z/R 

(3-13) 

9  9  9  2 

ill  m  z  +  n  +  e  i  -  L/R 

(3-14) 

The  derivation  in  Appendix  B  shows  that  wave  phase  is  given  by 

i  C.  -  q (9 , ip)  c.  -  q(60,O 

A(x,y)  -  A  (C,n)  -  Kpo»Rc2  J  c±  +  q(6,*)  Cj.  +  q(eQ,0  d? 

(3-15) 

where 

cx  -  [2(y-1)M2  +  4]/M2 

(3-16) 

c2  - (i  +  ^  *£>  Y~1 

(3-17) 

q(9,<))  -  1  +  (1-  3  cos26)/i|*3  +  1/4(1  +  3  cos29)/i(-6 

(3-18) 

29 


REPORT  MDC  A2582 
14  DECEMBER  1973 


GP74-03CXMS 


FIGURE  15 

COORDINATE  SYSTEM  FOR  MAIN-FLOW  SPHERE  PROBLEM 
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and 


cos  8  =  (c  cos  0Q  -  n  sin  6 Q) /t/j .  (3-19) 

2  2 

Equation  (3-15)  is  not  valid  if  at  any  point  in  space  V  ‘Mi  aa.  Since  the 
maximum  value  of  q  is  9/4  as  shown  in  Appendix  B,  we  must  have 

M2  <  J  or  M  <-|  (3-20) 

CO  9  00  3 

If  exceeds  this  limit,  V  will  reach  the  speed  of  sound,  and  a  shock  wave  will 
form. 


Figure  16  shows  plots  of  A' (0,n) /p®R  at  Mach  0.6  for  several  values  of  8q. 

The  plots  were  obtained  by  evaluating  Equation  (3-15)  numerically.  The  upper 
limit  used  was  l  ■  50,  which  effectively  is  the  same  as  l  m  <*>. 

These  calculations  of  A  are  based  on  a  compressible  but  nonvlscous  fluid. 

A  real  fluid  such  as  air  would  produce  a  flow  separation  aft  of  the  sphere  at 
very  low  Mach  numbers,  and  calculated  values  of  A  for  paths  traversing  the  flow 
separation  are  not  valid.  However,  the  calculations  do  provide  good  estimates  for 
forward-looking  paths.  In  fact  the  estimate  is  reasonable  for  forward-looking 
paths  even  if  M  >  2/3  since  the  shock  wave  will  occur  at  polar  angles  of  about  90° . 

Although  A  gives  a  complete  description  of  the  optical  degradation,  the  main- 
flow  0TF  may  also  be  required.  The  0TF  can  be  obtained  simply  by  applying 
Equation  (2-27)  to  A.  The  term  I  in  Equation  (2-27)  refers  to  the  beam  cross 
section.  Aircraft  are  streamlined,  and  therefore  do  not  have  as  severe  a  main- 
flow  optical  effect  as  the  sphere.  The  fact  that  the  sphere  does  produce  a  large 
effect,  however,  indicates  that  the  main-flow  effect  should  always  be  considered 
as  a  possible  source  of  optical  degradation. 

Shock  Waves 

In  this  subsection  a  very  brief  qualitative  description  of  the  optical  effect 
of  shock  waves  is  presented,  and  the  refraction  caused  by  a  plane  oblique  shock 
wave  is  computed. 

Qualitative  Description.  Shock  waves  are  very  thin  pressure  gradients  in 
aircraft  boundary  flow.  They  occur  only  when  some  part  of  the  flow  field  has  a 
speed  relative  to  the  aircraft  greater  than  the  speed  of  sound.  The  shock  wave 
thickness  is  only  a  small  fraction  of  a  centimeter,12  so  that  for  our  purposes  the 
shock  wave  may  be  considered  a  line  separating  two  regions  of  different  pressures, 
velocities,  and  densities.  Since  refractive  index  is  a  linear  function  of  density, 
the  shock  wave  is  also  the  interface  between  regions  of  different  refractive 
Indices.  Light  passing  through  a  shock  wave  will  be  refracted  according  to  Snell's 
law.  This  effect  is  generally  unimportant  in  aerial  photography,  since  the  image 
is  not  degraded  but  only  displaced  slightly.  An  unusual  exception  is  the  case  of 
the  Refaired  Nose  RF-4C  flying  at  high  subsonic  Mach  numbers.  As  air  moves  past 
the  nose  it  is  accelerated  to  a  speed  greater  than  the  speed  of  sound.  Under 
certain  conditions  a  shock  wave  form*  at  the  camera  window  as  shown  in  Figure  17. 
When  this  happens  a  fuzzy  band  appears  on  the  image  format  as  shown  in  the  figure. 
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FIGURE  16 

NORMALIZED  WAVE  PHASE  [A'{0,  rjJ/PooR]  vi  NORMALIZED  DISPLACEMENT  (r?) 

FOR  THE  MAIN-FLOW  SPHERE  PROBLEM 
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FIGURE  17 

FORMATION  OF  SHOCK  WAVE  ON  CAMERA  WINDOW  OF  THE  REFAIRED-NOSE 
RF-4C  AIRCRAFT  AND  ITS  EFFECT  ON  THE  IMAGING  SYSTEM 
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Light  from  the  ground  area  corresponding  to  the  fuzzy  band  arrives  at  the  camera 
with  an  angle  such  that  part  of  the  light  does  not  cross  the  shock  wave  but  the 
rest  does  cross  and  is  therefore  refracted.  The  result  is  a  double  image  which  is 
the  fuzzy  band. 

The  description  of  the  shock  wave  as  the  interface  between  media  of  different 
refractive  indices  is  somewhat  oversimplified.  In  general  shock  waves  are  curved, 
and  the  refractive  index  on  either  side  varies  continuously  along  the  shock  wave. 
This  complication  causes  the  shock  wave  to  not  only  refract  light  passing  through 
it,  but  also  to  distort  the  optical  wavefront.  No  attempt  is  made  here  to  analyze 
such  a  distortion.  The  angular  deviation  caused  by  refraction  for  a  plane  wave  is 
analyzed,  however,  in  the  next  subsection. 


Refraction  by  a,  Plane  Shock  Wave.  Although  refraction  is  usually  not  a 
problem  in  aerial  photography,  it  is  Important  in  optical  systems  for  which  the 
target  must  be  accurately  located.  Here  we  derive  the  angular  refraction  caused 
by  the  plane  shock  wave  which  is  produced  by  the  geometry  shown  in  Figure  18. 

Air  travels  at  some  supersonic  speed  u^  parallel  to  the  horizontal  surface.  When 
the  air  reaches  the  ramp  inclined  at  an  angle  6,  the  air  is  turned  parallel  to  the 
ramp  and  a  shock  wave  is  formed.  We  define  the  angle  of  the  shock  wave  to  be  £5 
and  the  horizontal  and  vertical  components  of  the  deflected  stream  of  air  to  be 
i>2  and  v£  respectively.  The  air  densities  to  the  left  and  to  the  right  of  the 
shock  wave  are  and  P2  respectively.  The  temperature  and  speed  of  sound  to  the 
left  are  T®  and  a«o  respectively. 

We  will  show  how  P2  can  be  found  from  a  knowledge  of  u.,  p-^,  T»,  and  either 
6  or  6-  Once  p^  and  P2  are  known  it  is  a  simple  matter  to  find  their  corresponding 
refractive  indices  n^  and  n2«  Snell’s  law  then  gives  the  angle  of  refraction  for 
light  crossing  the  shock  wave. 


Before  going  into  the  derivation  it  is  necessary  to  define  the  critical  speed 
of  sound  ac.  Going  back  to  the  reference  system  of  Figure  13,  ac  is  defined  as  the 
freestream  speed  U  when  U  is  equal  to  the  freestream  speed  of  sound  a,,.  In  other 
words,  when  the  freestream  Mach  nunfcer  is  unity  -  U/a,,,  ■  1),  then  we  have  ac  » 
U  ■  a*,.  The  speed  ac  may  be  expressed  in  terms  of  U  and  M»  byll 


2  Y  -  1  +  2/\ 
1  r  +  1 


(3-21) 


> 


Notice  that  when  M,,,  ■  1,  Equation  (3-21)  reduces  to  ac  ■  U  as  expected. 

Going  back  to  Figure  18,  we  see  that  simple  geometry  provides  the  relationship 


tan  0  ■  — -  • 
2 


(3-22) 


1 
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FIGURE  18 

FORMATION  OF  A  PLANE  OBLIQUE  SHOCK  WAVE 


REPORT  MDC  A2582 
14  DECEMBER  1973 


By  using  the  lews  of  conservation  of  mass,  momentum,  and  energy,  it  can  be  shown 
that  U2  and  V£  must  also  satisfy  the  following  equation. H 


v 


2 

2 


Equations  (3-22)  and 


(ul  "  “2)2  (ul“2  ~  ac> 
y+1  U1  ulu2  +  ac 
(3-23)  may  be  rewritten  as 

v'2 

tan  0  “  — r 


and 


where 


(“{  “  up2  (u£u£  -  1) 
r+T(u^)  -  V2  +  1 


(3-23) 


(3-24) 


(3-25) 


(3-26) 


Figure  19  plots  v^  versus  U£  for  both  Equations  (3-24)  and  (3-25).  The  values 
of  0  and  u£  for  this  figure  are  respectively  15*  and  2.  The  two  curves  intersect 
at  three  points  C,  D,  and  E.  Values  of  u£  and  v£  corresponding  to  these  points 
satisfy  both  equations.  The  value  of  u£  at  point  A  on  the  u^-axls  is  equal  to  u£. 
Therefore  values  of  (u£,  v£)  to  the  right  of  point  A  require  that  the  flow  Increase 
speed  as  it  passes  through  the  shock  wave.  Such  an  Increase  can  be  shown  to  be 
physically  impossible  since  entropy  would  be  required  to  decrease. H  We  may  therefore 
eliminate  the  values  of  u£  and  V2  at  point  E  from  our  set  of  solutions.  The 
remaining  two  points,  C  and  D,  represent  the  possible  velocity  components  for  the 
flow  to  the  right  of  the  shock  wave.  Point  C  is  said  to  represent  the  solution 
for  a  weak  shock  wave  and  D  a  strong  shock  wave,  since  the  velocity  change  through 
the  two  shock  waves  are  respectively  small  and  large.  Generally  the  shock  wave 
produced  experimentally  will  be  the  weak  shock  wave . H 

Figure  19  can  be  used  to  find  not  only  u£  and  v£  but  also  the  shock  wave 
angle  8.  Again  by  using  the  conservation  laws,  it  can  be  shown  that^ 


or 


tan  8 


(3-27) 


tan  8 


(3-28) 


FIGURE  19 

GRAPHICAL  SOLUTION  OF  NORMALIZED  VELOCITY  COMPONENTS 
(u'2.  v'2)  FOR  THE  PLANE-SHOCK-WAVE  PROBLEM 
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To  find  Che  angle  S  of  Che  weak  shock  wave,  a  verclcal  line  Is  drawn  Chrough 
polnC  C  as  shown  in  Figure  20.  This  line  Intersects  Che  U2~axis  aC  point  F. 

Since  Che  lengchs  of  Che  line  segmenCs  FA  and  CF  are  ui  -  u£  and  vo  respecCively , 
Che  angle  FCA  is  equal  Co  6.  Going  a  liCCle  fur Cher,  a  perpendicular  Co  Che 
line  defined  by  A  and  C  can  be  drawn  chrough  Che  origin  0.  The  angle  between  this 
perpendicular  and  the  u^-axis  is  also  equal  Co  0.  Figure  20  is  now  a  construction 
of  Che  actual  ramp  and  the  weak  shock  wave.  A  similar  construction  for  Che  strong 
shock  wave  would  show  Chat  Che  strong  shock  wave  has  a  larger  shock  wave  angle. 

In  Che  limit  as  Che  ramp  angle  8  approaches  zero,  Che  strong  shock  wave  angle 
approaches  90°.  On  Che  other  hand,  Che  weak  shock  wave  angle  approaches  a  value 
given  by 


sin  g 


_1_ 

M 

00 


(3-29) 


where  M,,,  ■  U/a^,  «  u^/a^.  Equation  (3-29)  can  be  derived  by  combining  Equations 
(3-21),  (3-23)  and  (3-27)  and  Caking  the  limit  u^uj.  The  value  of  a^  can  be 
found  by  using  Equation  (3-12). 

Figure  21  shows  chat  as  the  ramp  angle  8  increases,  points  C  and  D  approach 
each  ocher  until  they  coincide  at  Che  maximum  ramp  angle  emax*  lf  the  flow 
deflection  produced  by  an  obstacle  is  greater  than  8may ,  a  plane  shock  wave  cannot 
be  formed;  there  is  no  real  solution  to  the  set  of  equations,  (3-24)  and  (3-25). 
Instead  a  curved  detached  bow  wave  is  formed  as  shown  in  Figure  22. 

After  the  values  of  uj,  v£»  and  g  have  been  obtained,  the  density  change 
across  the  shock  wave  can  be  found  by  using  the  equation  for  mass  conservation, H 

PjU^  sin  g  -  p 2(u2  sin  ®  “  v2  COS  (3-30) 


or 


u^  sin  g 

P1  u'  sin  g  -  v'  cos  g 


(3-31) 


We  see  now  that  the  density  P2  can  be  obtained  by  finding  U2>  v2»  and  B 
using  Equations  (3-24)  and  (3-25)  and  then  by  applying  Equation  (3-31).  If  the 
angle  8  is  known  a  priori,  however,  p2  can  be  found  without  solving  for  U2  and  V2 
first.  As  shown  in  Appendix  C,  the  density  P2  may  also  be  expressed  as 


p2 


P1 


(y+1)M2  sin2  g 

0D 

(y-l)M2  sin2  g+2 
00 


(3-32) 


The  density  change  is 


2M2  sin2  8  -  2 

Ap  -  p,  -  p  -  p  - - - = -  ,  (3-33) 

1  (y-l)JT  sin  8  +  2 


last  commev 
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FIGURE  20 

GRAPHICAL  SOLUTION  OF  SHOCK-WAVE  ANGLE  0  FOR  THE  PLANE- 
SHOCK-WAVE  PROBLEM 
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FIGURE  21 

UPPER  LIMIT  OF  RAMP  ANGLE  0MAX  FOR  A  PLANE  SHOCK  WAVE 
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and  the  refractive  index  change  is 

An  -  n2  -  nx  -  (1  +  Kp^  -  (1  +  Kp2>  -  KAp  (3-34) 

or 

M2  sin2  6-1 

An  ■  2Kp.  — - - - - -  (3-35) 

(y-1)*T  sin  6  +  2 

where  K  is  given  by  Equation  (2-6)  or  approximately  by  K  ■  0.223  x  10  m  /kg. 

We  now  wish  to  find  the  refraction  angle  A<f  for  a  beam  intersecting  the  shock 
wave.  We  let  if  be  the  angle  between  the  beam  and  the  normal  to  the  shock  wave 
as  shown  in  Figure  23.  Snell's  law,  n  sin  if  *  const.,  is  differentiated, 

dn  sin  if  +  n  cos  if  d$  ■  0  (3-36) 

or 

dif  -  -  tan  if  dn.  (3-37) 

Since  An  is  very  small,  Equation  (3-37)  may  be  written  as  a  difference  equation, 

Aif  ■  -  tan  if  An.  (3-38) 

Finally,  Equation  (3-35)  is  substituted  in  Equation  (3-38). 

2  2 

fr  sin  6-1 

A<f  -  —  2Kp.  — - - - - -  tan  if  (3-39) 

(y-l)Mz  sin  6+2 
00 

The  minus  sign  indicates  that  when  the  beam  crosses  the  shock  wave  from  Region  I 
(lower  density)  to  Region  II  (higher  density),  the  beam  is  refracted  toward  the 
shock  wave  normal  as  shown  in  Figure  23.  Figure  24  shows  Aif/p^  tan  if  plotted  as 
a  function  of  6  for  various  Mach  numbers.  The  dotted  line  divides  the  plot  into 
strong  and  weak  shock  waves.  Figure  25  shows  A$/p.  tan  if  plotted  as  a  function  of 
M«,  for  both  the  strong  normal  shock  wave  (6  ■  0°,  6  *  90°)  and  for  the  shock  wave 
at  0  *  9max*  Since  strong  attached  shock  waves  generally  do  not  occur  in  practice, 
the  curve  for  9  *  8max  represents  the  maximum  angular  deviation  that  is  likely  to 
occur.  The  curve  for  9  M  0°  and  6  ■  90°  represents  the  maximum  angular  deviation 
that  can  occur  under  any  circumstance. 
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FIGURE  23 

REFRACTION  OF  A  BEAM  BY  A  PLANE  SHOCK  WAVE 
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FIGURE  25 

NORMALIZED  REFRACTION  ANGLE  (Mtpy  TAN  <fi)  vs  FREE-STREAM  MACH  NUMBER 
FOR  THE  STRONG  NORMAL  SHOCK  WAVE  AND  THE  SHOCK  WAVE  AT 
THE  RAMP-ANGLE  LIMIT  (0  -  0MAX) 
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IV.  Random  Effects 

The  purpose  of  this  section  is  to  present  imperical  methods  for  estimating  the 
optical  effects  of  turbulent  boundary  layers  and  separated  flow.  First  analytic 
equations  for  the  average  OTF  and  the  variance  of  the  wavefront  phase  a 2  are  pre¬ 
sented.  After  a  general  description  of  the  turbulent  boundary  layer,  the  Stine- 
Winovich  experiment  is  described.  The  results  of  the  experiment  are  uaed  to 
evaluate  o^  and  the  average  OTF  for  a  TBL.  Separated  flow  is  presented  in  a 
similar  way.  After  a  general  description  of  separated  flow,  data  from  the  RF-4C 
wind  tunnel  test  is  used  to  evaluate  o*  and  the  average  OTF.  Finally,  Kelsall's 
experiment  with  a  shearing  interferometer  is  described.  This  experiment  is  impor¬ 
tant  because  the  shearing  interferometer  technique  is  capable  of  measuring  MTF's 
directly. 

Equations  for  Average  OTF  and  a £ 


The  purpose  of  this  subsection  is  to  present  analytic  equations  for  the 
average  OTF  and  for  o^. 


Equation  for  Average  OTF.  Because  the  aperture  wave  phase  4  is  an  unknown 
random  variable,  the  OTF  is  also  an  unknown  random  variable.  It  is  possible,  how¬ 
ever,  to  derive  a  model  for  the  average  OTF  based  on  the  stochastic  properties  of 
A.  The  derivation  of  the  model  is  given  by  O'NeillS  and,  in  more  complete  form,  by 
Barakat.  3  Since  both  of  these  derivations  make  use  of  rather  advanced  theorems 
from  probability  theory,  a  third  simplified  derivation  is  given  in  Appendix  D.  All 
of  these  derivations  assume  that  A  defines  a  spatially  stationary  Gaussian  random 
process . 


As  stated  by  Barakat: 

"This  is  purely  a  working  hypothesis  backed  in  part  by  a  central  limit 
theorem  argument.  The  spatial  statlonarlty  condition  ...  is  mainly  one 
of  convenience  and  can  be  relaxed  at  the  cost  of  more  encumbering 
mathematics." 

The  derivations  show  that  the  expected  value  of  the  OTF  is 

E[T(ff‘Xf)]“  Tc(u’Xf  )  exp  {-k2[02-*<x,y)]}  (4-1) 

where  $(x,y)  and  are  respectively  the  covariance  function  and  variance  of  A. 

The  fraction  tc  is  the  OTF  which  would  result  if  there  were  controlled  aberrations 
only.  If  there  are  no  controlled  aberrations,  tq  becomes  the  diffraction  OTF,  and 
Equation  (4-1)  becomes 

E[T(xf*Xf)]  "  To(  Xf-'Xf-  )  exp{-kV-<Kx,y)]}  .  (4-2) 

Throughout  this  section,  only  random  aberrations  are  considered,  and  Equation  (4-2) 
is  used  Instead  of  Equation  (4-1). 
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If  we  define  the  ratio  of  t  to  t q  (or  t  q)  as  the  OTF  of  the  random  fluctuations 
only,  then  the  expected  value  of  the  random  OTF,  tr,  is 

E  [Tr(#-U  )  1  •  “P<  -K2[»2-«(x,y>!  }.  (4-3) 

By  definition  we  have 

*  -  E(A(x',y')A(x'+x,  y'+y) ]  -  E[A(x',y') ]E[A(x'+x,  y'+y)]  (4-4) 

and 

o2  -  E  [{A(x',y')  -  E[A(x',y')])2]  .  (4-5) 

Expanding  Equation  (4-5)  yields 


a2  -  E[A2(x',y')]  -  (E[A(x',y')]}2 


*(0,0). 


(4-6) 


Because  spatial  stationarity  is  assumed,  $  is  a  function  of  x  and  y  only,  E(A)  is 
constant,  and  a 2  is  constant.  Without  loss  of  generality  we  take  the  expected 
value  of  A  to  be  zero.  We  now  have 


i 


l 


*(x,y)  -  E[A(x',y')A (x'+x,y'+y) ] 

(4-7) 

and 

o2  -  E[A2(x',y')]  -  *(0,0). 

(4-8) 

Since  obtaining  averages  from  an  ensenfcle  is  impractical,  time 
under  the  assumption  that  time  and  ensemble  averages  are  equal 
(4-3),  (4-4),  and  (4-8)  become 

averages  are  used 
.  Equations  (4-2) , 

(4-9) 

•*p{-k2I»2-»Cx.y)]} 

(4-10) 

*(x,y)  -  A(x',y')A(x'+x,y"+y) 

(4-11) 

and 

o2  ■  *(0,0)  ■  A2(x',y') 

(4-12) 

where  the  ovetbar  indicates  time  average.  Since  and  Tq  are 

both  real,  and  T 

are  average  MTF's  as  well  as  average  OTF's. 


It  is  still  necessary  to  find  a  functional  form  for  * .  Again  we  adopt  the 
attitude  of  Barakat.  (Barakat's  Rz  is  equivalent  to  *) . 
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"...in  lieu  of  any  contrary  evidence  Rz  is  taken  to  be  Gaussian  as  a 
mathematical  convenience  rather  than  because  it  has  a  physical  founda¬ 
tion.  I  have  reason  to  believe  that  the  actual  shape  of  Rz  (provided 
it  be  a  smooth  function)  is  less  important  than  its  spatial  spread." 


The  function  $(x,y)  can  now  be  written 


4>(x,y)  -  a2  expj^-^r  + 


U- 


3] 


(4-13) 


where  ax  and  ay  are  the  correlation  lengths  in  the  x  and  y  directions  respectively. 
Throughout  the  rest  of  this  report,  the  random  wavefront  aberrations  are  assumed  to 
be  isotropic,  that  is. 


Equation  (3-12)  may  now  be  written 

<K(x,y) 

The  average  random  OTF  becomes 


4>(r) 


2  -(r/a)‘ 
o  e 


(u)  "  exp{  *k2a2  [l-e'(r/a)2]  } 


or 


where 


TrM 


exp  < 


r 

Xf’ 


(4-14) 


(4-15) 


(4-16) 


(4-17) 


(4-18) 


We  now  have  a  two  parameter  (a2  and  a)  model  for  the  time  average  OTF.  Notice  that 
Tr  never  falls  to  zero  but  has  a  minimum  value  given  by 


.2  2 

-k  a 


min 


(4-19) 


Equation  (4-16)  is  plotted  in  Figure  26  for  different  values  of  o/X  and  a.  The  a 
parameter  determines  how  low  Tr  will  drop,  and  the  correlation  length  a  is  a  scale 
factor  for  the  spatial  frequency  v.  Usually  a2  and  a  must  be  found  experimentally; 
however,  the  following  analytic  development  does  show  the  dependence  of  o2  on  such 
parameters  as  density  and  velocity. 
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fp  for  Different  Values  of  ofk 
a  is  Constant 


rR  for  Different  Values  of  a 
oA  =  0.2 


FIGURE  26 

AVERAGE  RANDOM  OTF,  rR,  FOR  VARIOUS  VALUES  OF  o/X  AND  a 
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Equation  for  o£.  Both  Chernov1^  and  Tatarskii  have  derived  an  equation  for 
02  in  terms  of  the  variance  of  the  refractive  index  as  a  special  case  of  a  more 
general  problem.  Chernov  has  also  used  a  simple  and  direct  epproach,  but  has 
omitted  some  of  the  steps.  For  this  reason,  a  complete  direct  derivation  is  given 
In  Appendix  E.  The  equation  for  o2  is  found  to  be 


a2  ■ 


(4-20) 


where  a  is  wavefront  correlation  length,  L  is  optical  path  length  through  the 
medium,  and  u?  is  the  variance  of  the  refractive  index.  The  definition  of  u  is 


n  -  n. 


so  that 


^  ■  (n-n)2 


n2-n2. 


(4-21) 


(4-22) 


The  assumptions  made  in  the  derivation  are  that  y<<l  and  that  y  has  a  Gaussian  co- 
variance  function. 


y  (x',y '■,z')u(x>+x,y''+y,z'+z) 


U1  exp[-(x2+y*>z2)/a2] 


(4-23) 


where  a<<L.  The  parameter  a  in  Equation  (4-23)  is  the  same  as  the  wavefront  corre¬ 
lation  length.  In  fact,  it  is  shown  in  Appendix  E,  that  Equation  (4-15),  which 
gives  the  functional  form  for  9,  is  implied  by  Equation  (4-23). 

One  way  of  estimating  y2  is  to  relate  the  refractive  index  fluctuation  to 
density  fluctuations  and  then  the  density  fluctuations  to  pressure  fluctuations, 
which  can  be  measured  at  the  airplane  skin.  The  transfer  from  density  to  pressure 
fluctuation  is  made  by  assuming  that  the  pressure  changes  are  adiabatic.  We  start 
by  substituting  Equation  (2-5)  into  Equation  (4-22) 


y2  ■  n2-  Y&  ■  (1  +  Kp)2  -  1+Kp 
y2  -  1  +  2Kp  +  K2p2  -  (1  +  2Kp  +  K2p2) 
y2  -  K2(^  -  p2)  -  K2 (p  -  p) 1  -  K2(Ao)2 


(4-24) 


(4-25) 


(4-26) 


where  p  is  density  and  K  is  a  constant  given  by  Equation  (2-6)  or  approximately  by 
K  ■  0.223  x  10* 3  m3/kg.  The  term  Ap  represents  density  variation  from  the  mean. 
The  adiabatic  law, 


Pp"Y  - 


const 


(4-27) 


where  P  is  pressure  and  y  is  the  ratio  of  specific  heats  (-1.4),  is  differentiated 
to  give 


p"YdP  -  YPp'<Y+1)dp  -  0. 
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Since  we  are  working  with  small  changes  in  pressure  and  density,  we  can  write  this 
differential  equation  as  a  difference  equation.  Si&stituting  AP  and  Ap  for  dP  and 
dp  respectively  and  rearranging.  Equation  (4-28)  becomes 


Ap  “  -%■  AP. 
rP 


(4-29) 


Generally  pressure  fluctuation  data  are  given  in  terms  of  the  fluctuating  pressure 
coefficient,  which  is  defined  by 


/  AP2 


Cp  |p.u2 


(4-30) 


where  U  is  the  airplane  speed  and  P®  is  the  freestream  density.  Combining 
Equations  (4-20),  (4-26),  (4-29),  and  (4-30)  we  have 


2  /T  K2  .  P.4U4  _  2 

a  "  7  7  "*•  “P“  ®P 


(4-31) 


where  P®  is  the  free- stream  pressure.  This  equation  can  be  simplified  by  using  mach 
number  Instead  of  speed. 


H»  -  0/s*. 

The  speed  of  sound  a.  is  given  by  Equation  (3-12) 


/yET 

"  /  P. 


^287yT„ 


(4-32) 


(4-33) 


where  T®  is  freestream  temperature  and  all  units  are  MRS. 
Substituting  Equations  (4-32)  and  (4-33)  into  Equation  (4-31)  yields 


V2n  2  -  „  4  2 

^  Poo  fiLMoO  Cp  s 


(4-34) 


Equation  (4-34)  predicts  the  wavefront  variance  when  the  fluctuating  pressure  coef¬ 
ficient  Cp  is  known.  It  will  be  shown  shortly  how  this  equation  can  be  applied  to 
the  TBL  problem. 

Turbulent  Boundary  Layer 

The  purpose  of  this  subsection  is  to  present  a  brief  description  of  the  turbu¬ 
lent  boundary  layer  and  then  to  apply  the  Stine-Wlnovlch  experimental  data  to  the 
average  OTP  and  o?  models  (Equations  (4-17)  and  (4-34)). 

Description.  Consider  a  thin  flat  plat  immersed  in  a  fluid  moving  parallel  to 
the  plate  surface  at  speed  U  as  shown  in  Figure  27.  A  Cartesian  coordinate  system 
is  defined  by  letting  the  x-axis  lie  along  the  top  of  the  plate  in  the  direction  of 
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LBL  =  Laminar  Boundary  Layer 
TBL  •  Turbulent  Boundary  Layer 


FIGURE  27 

TURBULENT  BOUNDARY  LAYER  ON  A  FLAT  PLATE 
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Che  fluid  flow  and  by  letting  the  y-axis  be  perpendicular  to  the  plate,  intersecting 
the  x-axis  at  the  leading  edge.  We  let  u  and  v  represent  the  x  and  y  fluid  velocity 
components  respectively.  From  the  definition  of  the  problem  the  velocity  components 
have  the  values 


u  ■  U  and  v  ■  0  (4-35) 

at  large  distances  from  the  plate.  On  the  surface  of  the  plate,  however,  friction 
imposes  the  "no  slip"  condition  on  the  fluid, 

u  -  v  *  0.  (4-36) 


Most  of  the  necessary  velocity  change  takes  place  in  a  very  narrow  region  just  out¬ 
side  the  plate  surface.  This  region  is  called  the  boundary  layer.  In  Figure  27  the 
thickness  of  the  boundary  layer  is  greatly  exaggerated  and  the  boundary  layer  on  the 
bottom  of  the  plate  is  not  shown.  Near  the  leading  edge  of  the  plate,  the  velocity 
change  is  smooth  with  distance  away  from  the  plate,  and  the  vertical  velocity  com¬ 
ponent  is  negligible.  Under  these  conditions  we  say  that  the  boundary  layer  is 
laminar.  This  laminar  boundary  layer  gradually  increases  In  thickness  with  distance 
from  the  leading  edge.  At  some  point  behind  the  leading  edge,  the  boundary  layer 
may  become  turbulent.  Random  variations  occur  in  the  velocity  components  u  and  v 
and  in  the  density  p.  In  addition  the  rate  of  thickness  increase  with  distance  from 
the  leading  edge  is  larger  than  that  of  the  laminar  boundary  layer. 


The  transition  from  laminar  to  turbulent  occurs  at  a  position  corresponding  to 
a  Reynolds  nunber  in  the  range 

3  x  105  <  R  <  3  x  106  (4-37) 

x 

where  the  characteristic  length  associated  with  the  Reynolds  number  is  the  distance 
from  the  leading  edge.*0  All  other  parameters  associated  with  the  Reynolds  number 
are  taken  at  freestream.  The  Reynolds  number  is  given  by 


R 


x 


P»UX 

Poo 


(4-38) 


where  X  is  distance  from  the  leading  edge  and  Poo  is  the  freestream  coefficient  of 
viscosity.  Since  transonic  aircraft  have  Reynolds  numbers  on  the  order  of  1  x  10° 
per  meter,  the  boundary  layer  will  be  turbulent  almost  everywhere. 


The  effect  of  friction  on  the  speed  of  the  air  is  greatest  at  the  interface 
between  the  plate  and  the  air,  and  this  effect  decreases  rapidly  with  distance  from 
the  plate.  The  question  arises  as  to  hew  to  define  the  boundary  layer  thickness  6. 
The  general  practice  is  to  define  6  as  that  distance  from  the  plate  at  which  the 
average  horizontal  velocity  component  is  equal  to  99%  of  the  main-stream  speed. 

That  is. 


u (y  -  5)  -  0.99V. 


(4-39) 
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Although  £  Is  called  che  boundary  layer  thickness,  other  types  of  thicknesses  have 
been  defined.  One  example  is  the  displacement  thickness  £*,  which  is  defined  by 


6* 


c  V 

-J  1_^U  dy 


(4-40) 


where  Py  is  the  average  density  inside  the  boundary  layer  at  the  distance  y  above 
the  plate.  Falkner1**  has  proposed  a  model  for  turbulent  boundary  layers  with 


Reynolds  numbers  between  10 ^  and 
respectively  by 


10 


10 


According  to  this  model,  £  and  £*  are  given 


£  - 


0.1285  X 


*x 


1/7 


and 


6* 


0.0214  X  1  „ 

'  0  1/7  "  6  5  * 


(4-41) 


(4-42) 


Equation  (4-41)  enables  us  to  estimate  the  optical  path  length  L  through  the 
TBL  by  setting  it  equal  to  £. 


L  ■  5 


(4-43) 


Likewise,  Equation  (4-42)  enables  us  to  estimate  the  correlation  length  a. 
According  to  Lyon17  the  spatial  scale  of  the  fluctuating  pressure  on  the  skin  of 
aircraft  is  approximately  equal  to  the  boundary  layer  displacement  thickness  £*. 
Since  the  density  variations  are  assumed  to  be  directly  related  to  the  pressure 
variations  according  to  Equation  (4-29) ,  the  fluctuating~density  spatial  scale  and 
the  fluctuatlng-ref ractive-lndex  spatial  scale  are  also  approximately  equal  to  £*. 
Figure  28,  which  plots  the  Gaussian  covariance  function  given  by  Equation  (4-15), 
shows  that  the  spatial  scale  is  approximately  equal  to  3a;  that  is, 


or 


£*  -  3a 


i 

18' 


(4-44) 


(4-45) 


Now  from  Equations  (2-54),  (4-38),  (4-41),  (4-43),  and  (4-45)  we  can  write  the  fol¬ 
lowing  equations  for  L  and  a. 


L  -  0.0147  X6/7T„3/28pa,~1/7lf1/7 


a  ■  L/18 


(4-46) 

(4-47) 


All  units  are  MRS. 
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FIGURE  28 

GAUSSIAN  AUTOCORELATION  FUNCTION 
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It  must  be  remembered  here  that  Equation  (4-46)  applies  only  to  flat  plates, 
and  even  so  it  is  only  an  approximation.  Since  analyses  on  other  shapes  are  not 
available,  however,  this  equation  must  serve  as  an  estimate  for  aircraft  when  no 
other  information  about  the  TBL  thickness  is  available.  On  shapes  other  than  a 
flat  plate  the  distance  X  would  be  the  distance  along  the  skin  from  the  leading 
edge  or  point. 

Although  Equations  (4-46)  and  (4-47)  give  estimates  of  L  and  a.  Equation  (4-34) 
shows  that  Cp  is  still  needed.  Fortunately  Cp  is  known  to  be  fairly  constant  for 
subsonic  aircraft  TBL's.  The  value  quoted  by  Lyon^  and  confirmed  by  F-4  flight 
tests  is 

Cp  -  0.006.  (4-48) 

Apparently  we  now  have  a  method  of  evaluating  and  a  by  using  Equations  (4-34) , 
(4-46),  (4-47),  and  (4-48).  The  results  of  the  Stine-Winovich  experiment,  which 
will  be  described  presently,  show  that  this  method  is  essentially  correct  but  that 
some  changes  in  the  equations  are  necessary.  Readers  interested  in  a  more  detailed 
description  of  turbulent  boundary  layers  are  referred  to  Reference  10. 

Stine-Winovich  Experiment.5  The  purpose  of  this  wind  tunnel  experiment  was  to 
measure  the  optical  effect  of  the  turbulent  boundary  layers  which  formed  on  the  wind 
tunnel  walls.  The  experiment  was  conducted  at  the  Ames  Aeronautical  Laboratory, 
Moffet  Field,  California  in  the  one~by-three-foot  supersonic  tunnel  No.  1.  Figure 
29  shows  the  experimental  arrangement.  Light  from  a  tungsten-filament  source  is 
collected  by  a  lens  and  passes  through  a  0.00025-inch  diameter  source  aperture  at 
the  focus  point.  The  light  is  collimated  by  an  f/8.6  Newtonian  telescope  with  a 
21.45  inch  focal  length  and  projected  across  the  wind  tunnel  through  plane  windows. 
On  the  other  side  of  the  tunnel  the  light  is  focused  by  a  second  Newtonian  tele¬ 
scope,  identical  to  the  first.  One  of  several  circular  collection  apertures  is 
placed  at  the  focus  point  of  the  second  telescope.  The  radiant  power  of  the  light 
passing  through  the  collection  aperture  is  measured  with  a  photometer  and  the 
intensity  of  the  tungsten  source  is  monitored  with  an  electronic  photometer.  The 
TBL  thickness  was  found  by  using  a  hot  wire  anemometer.  Measurements  were  made  for 
the  normal  TBL  with  a  thickness  of  about  1/10  foot  and  for  an  artificially  thickened 
TBL  with  a  thickness  of  about  2/10  foot.  The  thickening  was  accomplished  by  apply¬ 
ing  a  roughness  coating  to  the  side  walls  of  the  tunnel  to  within  40  inches  of  the 
windows.  Figure  30  shows  the  positions  of  the  anemometer  probe  and  the  roughness 
coating  relative  to  the  window. 

The  analytic  approach  taken  by  Stine  and  Winovich  was,  in  this  writer's 
opinion,  inadequate.  However,  there  is  no  reason  why  their  experimental  data 
should  not  be  used  to  examine  our  model  for  the  average  OTF  of  a  TBL. 


Appendix  F  shows  that  the  average  theoretical  ratio  Q  of  the  radiant  power 
emerging  from  the  collection  aperture  to  that  emerging  from  the  source  aperture  is 
given  by 


u  — 
rb 


(4-49) 
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FIGURE  29 


STINE-WINOVICH  EXPERIMENTAL  ARRANGEMENT 
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FIGURE  30 

VIEW  OF  WIND-TUNNEL  WALL  FOR  STINE-WINOVICH  EXPERIMENT 
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where 

D  ■  telescope  aperture  diameter 

£  -  telescope  focal  length 
ra  ■  radius  of  source  aperture 
rb  *  radius  of  collection  aperture 

=  first  order  Bessel  function  of  first  kind 

The  appendix  shows  that  Tr  is  squared  because  the  light  beam  must  traverse  two 
boundary  layers. 

The  problem  here  is  to  match  values  of  Q  computed  from  Equation  (4-49)  and  our 
model  for  Tr  with  the  values  of  Q  measured  in  the  experiment.  Values  of  Q  were 
measured  for  various  mach  numbers  and  densities  for  both  the  natural  boundary  layer 
and  the  artificially  thickened  boundary  layer.  Equation  (4-49)  was  evaluated  by 
using  Equations  (2-7),  (2-35),  (4-17),  (4-34),  (4-43),  (4-47),  and  (4-48).  Although 
the  Newtonian  telescopes  had  a  central  obscuration,  the  diameter  of  which  was  one- 
tenth  the  aperture  diameter,  the  diffraction  OTF  of  a  circle  can  be  used  with  negli¬ 
gible  error.  Figure  31  shows  a  comparison  of  the  diffraction  OTF's  for  a  circular 
aperture  and  for  the  center-blocked  circular  aperture. 

Figure  32  shows  both  the  measured  and  computed  values  of  Q  vs r^  for  the 
highest  density  and  mach  number.  The  computed  values  are  represented  by  the  dashed 
line.  Obviously  the  experimental  values  are  much  lower  than  the  corresponding  com¬ 
puted  values.  The  implication  is  that  the  actual Tr  is  much  lower  than  we  would 
predict.  A  good  match  can  be  achieved,  however,  with  only  minor  changes  to  our 
model.  Specifically,  the  coefficient  in  Equation  (3-47)  is  changed  from  1/18 
(0.056)  to  0.08  and  the  value  of  Cp  is  changed  from  0.006  to  0.05. 

a  -  0.086  -  0.08L  (4-50) 

Cp  -  0.05  (4-51) 

The  solid  line  shows  a  plot  of  vs.  r^  which  was  computed  by  using  Equations  (4-50) 
and  (4-51)  instead  of  Equations  (4-47)  and  (4-48).  This  plot  shows  good  agreement 
between  the  computed  and  measured  values  of  Q.  Figure  33  shows  plots  of  Q  vs.  rfc 
using  the  new  Equations  for  various  values  of  mach  number,  density  and  boundary 
layer  thickness. 

By  substituting  Equations  (2-7),  (4-43),  (4-50),  and  (4-51)  into  Equation 
(4-34) ,  the  variance  may  now  be  expressed  as 

o2  -  4.4  x  10”1262p0O2Ho4.  (4-52) 

All  units  are  MKS.  For  a  flat  plate  Equations  (4-32),  (4-33),  (4-43),  and  (4-46) 
may  be  substituted  into  Equation  (4-52)  to  yield 
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Normalized  Spatial  Frequency  r/D 


m 
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The  correlation  length  from  Equation  (4-50)  becomes 


a 


7.66  x 


10~4p, 


-l/7_  1/28  — 1/ 7„6 / 7 
„  I*  m,,,  a  . 


(4-54) 


Again  all  units  are  MKS . 

Even  though  the  value  of  Cp  had  to  be  changed  by  an  order  of  magnitude  to  make 
our  model  work,  the  equations  for  t j  (Equation  (4-17))  and  for  a  2  (Equation  (4-34)) 
are  apparently  reliable.  The  original  failure  of  the  model  should  not  be  judged 
too  harshly.  It  must  be  remembered  that  we  were  trying  to  describe  density  fluctua¬ 
tions  inside  the  wind  tunnel  boundary  layer  from  pressure  fluctuations  on  the  sur¬ 
face  of  an  aircraft. 

_  A  valid  criticism  of  the  Stine-Winovich  experiment  is  that  the  measured  value 
of  Q  may  be  affected  by  vibration.  The  fact  that  the  experimental  data  does  agree 
with  the  model  for  changes  in  mach  number,  density,  and  boundary  layer  thickness 
indicates  that  the  effect  of  vibration  may  well  be  negligible.  Even  if  vibration 
did  effect  the  experiment,  the  model  can  still  be  used  with  the  understanding  that 
0.05  is  an  upper  limit  on  Cp.  Besides  the  vibration  problem,  there  is  the  fact  that 
the  OTF  is  valid  only  for  monochromatic  light.  It  is  this  writer's  belief,  however, 
that  the  model  for  the  turbulent  boundary  layer  OTF  presented  in  this  section  is  the 
best  one  available.  An  accurate  model  must  await  direct  OTF  measurements  under  con¬ 
trolled  test  conditions. 

Separated  Flow 

Treatment  of  separated  flow  in  this  subsection  is  similar  to  the  treatment 
given  TBL's  in  the  last  subsection.  After  a  general  description  of  separated  flow, 
data  from  the  RF-4C  wind  tunnel  test  is  used  to  evaluate  and  a  for  the  average 
separated  flow  OTF. 

Description.  ^  Consider  the  surface  in  Figure  34  which  has  a  steady  convex 
curvature,  the  surface  curving  away  from  the  direction  of  flow.  In  this  situation 
the  mainstream  flow  is  retarded  and  the  pressure  increases  in  the  direction  of  the 
flow.  It  can  be  shown  that  static  pressure  is  essentially  constant  across  the 
boundary  layer;  that  is 


_3P 

3y 


0 


(4-55) 


inside  the  boundary  layer.  This  condition  implies  that  static  pressure  inside  the 
boundary  layer  increases  in  the  direction  of  flow. 


3P 

3x 


> 


0 


(4-56) 


The  fluid_element  ABCD  will  have  an  average  pressure  P  on  the  AD  side  and  an  average 
pressure  P  +  (3P/3x)  6x  on  the  BC  side.  The  resultant  pressure  force,  in  addition 
to  friction,  will  slow  the  element.  The  deceleration  is  greatest  near  the  surface 
where  the  element  is  remote  from  the  mainstream.  The  figure  shows  the  change  in  the 
average  velocity  profiles.  The  deceleration  of  the  fluid  implies  that  the  boundary 
layer  must  thicken  to  conserve  mass.  At  some  point  S  on  the  surface  the  velocity 
gradient  3u/5y  vanishes.  Beyond  this  point  reverse  flow  may  occur  as  shown  in  the 
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figure.  The  dotted  line  extending  from  S  in  the  downstream  direction  is  drawn  such 
that  the  mass  flow  rate  between  it  and  the  edge  of  the  boundary  layer  is  equal  to 
the  boundary  layer  mass  flow  rate  upstream  of  S.  The  boundary  layer  actually  sepa¬ 
rates  from  the  surface  with  the  region  below  the  dotted  line  being  the  separated- 
flow  region. 

RF-4C  Wind  Tunnel  Test.  The  purpose  of  this  test  was  to  measure  the  effect  of 
the  separated  flow  field  behind  the  RF-4C  nose  chin  on  photographic  resolution.  The 
test  was  conducted  in  the  transonic  Propulsion  Wind  Tunnel  16T  at  the  Arnold  Engi¬ 
neering  Development  Center,  Arnold  Air  Force  Station,  near  Tullahoma,  Tenn.  Figure 
35(a)  shows  the  experimental  setup.  A  model  of  the  bottom  half  of  the  RF-4C  nose 
with  a  Chicago  Aerial  CA-120  camera,  is  mounted  on  a  splitter  plate.  Both  full 
scale  and  fifth  scale  models  were  used.  The  camera  took  pictures  of  a  standard 
Air  Force  resolution  bar  target  through  a  collimator  mounted  in  the  wind  tunnel 
ceiling.  The  camera  was  equipped  with  a  Xenotar  6-inch  lens  and  the  film  used  was 
Kodak  3414.  The  target  contrast  was  1000:1.  The  external  square  aperture  stop  was 
mounted  diagonally  with  respect  to  the  direction  of  flow  and  with  respect  to  the 
resolution  bars  as  shown  in  Figure  35(b).  Runs  were  made  for  various  exposure 
times,  aperture  settings,  air  densities,  and  mach  numbers  to  determine  optimum 
exposure  times  and  aperture  settings  for  different  flight  conditions. 

Although  the  RF-4C  Wind  Tunnel  Test  was  not  designed  for  MTF  computations,  a 
portion  of  the  image  resolution  data  can  be  used  to  develop  simple  models  for  a?  and 
a  to  be  used  in  Equation  (4-17) . 

1.  Background.  In  Section  II  it  was  shown  that  the  image  of  a  sinusoidal  tar¬ 
get  of  modulation  Cq  is  a  sinusoid  of  modulation  C^  given  by 

C1  -  C0|t(v)|,0  <.  C0  <_  1  (4-57) 

where  | t (v) |  is  the  MTF  of  the  imaging  system.  Normally  the  system  would  include 
the  boundary  flow,  aperture,  and  lens.  In  this  particular  case,  however,  lens 
degradation  is  insignificant.  We  now  define  system  resolution  as  the  largest 
spatial  frequency  umax  at  which  the  image  can  be  resolved  for  the  given  value  of 
Co>  Becuase  the  film  in  the  camera  and  the  eye  of  the  individual  recording  the 
resolution  are  not  optically  perfect,  vmax  is  always  less  than  \>  ,  the  spatial 
frequency  at  which  T  vanishes.  0 

T (vQ)  -  0  (4-58) 


In  order  to  find  vMX  threshold  modulation  (TM)  is  introduced.  The  TM  is  defined 
as  the  lowest  value  of  Ci  for  which  a  target  of  a  given  spatial  frequency  can  be 
resolved.  Figure  36  shows  a  typical  TM  plotted  as  a  function  of  spatial  frequency. 
The  TM  is  a  function  of  the  film  and  the  human  recorder.  From  the  definition  of  TM 
and  Equation  (4-57)  we  see  that  as  long  as  the  product  Cq |t  (v) |  is  greater  than 
TM(v),  the  image  can  be  resolved.  At  v  -  vmax  we  have 


(4-60) 
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FIGURE  35 

EXPERIMENTAL  ARRANGEMENT  FOR  RF4C  WIND-TUNNEL  TEST 
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<  I 


If  the  MTF  is  factored  into  its  diffraction  and  random  parts ,  Equation  (4-60)  may 
be  written 


TR^vmax^ 


TM(v  ) 

max 

0  0  max 


(4-61) 


In  analyzing  the  RF-4C  test,  tq  and  TM  are  known  as  are  values  of  vmPy  for  some 
particular  conditions.  For  all  practical  purposes  the  target  modulation  for  the 
1000:1  target  is  one. 


CO 


1000  -  1  _ 
1000  +  1 


(4-62) 


Using  this  information  we  model  a2  and  a  in  Equation  (4-17)  for  tr  such  that 
Equation  (4-61)  is  satisfied.  Although  the  target  in  the  test  was  a  rectangular-bar 
target  instead  of  a  sinusoidal  target,  it  was  felt  that  the  accuracy  in  determining 
vmax  would  be  sufficient  for  this  analysis. 

Equation  (4-61)  shows  that  expressions  for  TM  and  t~  are  needed.  Lauroesch 
et  al.  have  found  that  a  good  model  for  the  TM  is 


TM(v)  -  b0  +  bjV2 


(4-63) 


where  bo  and  b^  are  constants.*®  The  bg  term  Is  not  a  function  of  the  film  but  of 
the  human  observer.  As  spatial  frequency  increases  from  zero,  the  film  MTF  and 
.  grain  noise  become  increasingly  important  as  indicated  by  the  b^v2  term.  The  value 
of  bQ  is  approximately  0.03.  The  constant  b^  for  Kodak  3414  film  is  found  by 
solving  Equation  (4-63)  for  bj  and  using  the  data  point  TM  »  0.231  for  v  »  250  ran" 
or  lines /mm. I® 


b±  -  [TM(v)  -  bQ]/v2  -  3.22  x  10”6mm2 


Equation  (4-63)  is  now 


TM(v)  -  0.03  +  3.22  x  10_6v2 


(4-64) 


(4-65) 


where  v  is  in  lines/mn.  The  expression  for  tq  is  found  by  using  Equations  (2-32) 
and  (2-25). 


(v  Xf  \  /  v  Xf\  D 


D 

JL 

Xf 


This  equation  can  be  simplified  by  recognizing  that  Dx  »  Dy  and  defining 


„  Xf  Xf 

B  “  —  ■  D”‘ 
x  y 


(4-66) 


(4-67) 
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Equation  (4-66)  becomes 


T0(vx,vy)  “  (1  *  Bvx5  (1  "  Bvy)f  0  —  Vx,Vy  -B 


(4-68) 


With  the  aperture  diagonal  to  the  target,  the  resolution  is  not  a  value  of 
either  vx  or  vy  but  of  v,  the  spatial  frequency  along  the  diagonal  between  the 
x  and  y  directions.  Since  displacement  along  the  diagonal  is  given  by 


r  -Vx+y  -  i/5"  x  ■  /2  y, 


(4-69) 


the  spatial  frequency  v  is  given  by 


v  -■Jv2  +  v2 
lx  y 


vT  v  -  /I  v 
x  y 


Substituting  Equation  (4-70)  into  Equation  (4-68)  we  have 


(4-70) 


t  (v)  -  1  -  /I  Bv  +  B2v2/2,  0  <  v  < 


(4-71) 


The  camera  aperture  setting  was  i/22, 


rr~  m  TT  m  22. 

D  D 
x  y 

Since  the  source  was  white  light,  the  wavelength  is  taken  as 
X  ■  0.55  x  10  2  mm  . 


(4-72) 


(4-73) 


Substituting  Equations  (4-72)  and  (4-73)  into  Equation  (4-67)  gives  B, 
B  ■  0.0121  mm, 
and  Equation  (4-71)  becomes 

tq (v)  -  1  -  0.0171  v  +  7.32  x  10~5v2. 


(4-74) 


(4-75) 


Equations  (4-65)  and  (4-75)  are  the  expressions  needed  for  TM  and  tq  respectively. 
Plots  of  these  equations  are  shown  in  Figure  37.  We  are  now  in  a  position  to 
look  at  the  test  data. 

2.  Data.  Most  of  the  data  collected  during  the  test  is  eliminated  from  the 
analysis  because  of  the  two  following  considerations.  First,  we  are  interested 
in  the  average  of  the  random  MTF,  and  therefore  only  data  related  to  the  longest 
exposure  times  for  each  configuration  is  considered.  Second,  it  is  necessary  to 
know  the  MTF  of  the  lens  if  the  lens  introduces  significant  degradation.  At  the 
i/22  aperture  setting,  performance  was  nearly  diffraction  limited,  but  at  smaller 
f  numbers  the  lens  produced  significant  degradation.  Since  the  MTF  of  the  lens 
at  the  various  f  numbers  is  unknown,  only  data  from  the  i/22  aperture  settings 
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are  analyzed.  Table  I  below  lists  the  resolution  data  considered. 1  Each  resolution 
given  Is  the  average  resolution  for  ten  exposures. 


Table  I  IMAGE  RESOLUTION 

(LINES /MM) 

FOR  RF-4C  WIND 

TUNNEL  TEST 

Scale 

Full 

Full 

1/5 

Pressure  altitude  (ft) 

8,000 

40,000 

8,000 

Exposure  time  (sec) 

0.010 

0.010 

0.001 

Mach  Number 

.7 

60 

56 

.8 

37 

.85 

39 

36 

.9 

40 

66 

36 

.95 

35 

47 

The  fifth-scale  40,000  ft  data  was  not  used  because  the  random  flow  degradation 
was  not  significant.  The  data  In  Table  1  shows  that  no  strong  correlation  exists 
between  Mach  number  and  resolution  for  Mach  numbers  between  .8  and  .95.  Since 
optical  degradation  is  greatest  in  this  Mach  range,  only  data  from  this  range  was 
used,  and  Mach  number  was  excluded  as  a  variable. 

A  possible  explanation  for  the  lack  of  correlation  may  be  that  there  was 
little  change  in  the  actual  air  speed  in  the  vicinity  of  the  flow  separation.  The 
shock  wave  in  front  of  the  flow  separation  shown  in  Figure  1  indicates  that  the 
air  has  been  accelerated  to  some  supersonic  speed  before  reaching  the  shock  wave. 

A  property  of  shock  waves  is  that  the  flow  deceleration  through  a  shock  wave 
increases  with  an  increase  of  flow  speed  into  the  shock  wave.  It  is  therefore 
reasonable  to  expect  that  changes  in  the  actual  flow  Mach  number  in  the  vicinity  of 
separated  flow  are  not  as  great  as  the  changes  in  the  freestream  Mach  number. 

Table  II  shows  the  data  used  in  the  analysis. 


Table  II  DATA  USED 

IN  ANALYSIS 

OF  RF-4C  WIND 

TUNNEL  TEST 

Scale 

Full 

Full 

1/5 

Pressure  Altitude 

8,000 

40,000 

8,000 

v  (mm  ^) 

max 

39 

66 

36 

0.444 

0.190 

0.479 

TM(v  ) 
max 

0.0305 

0.0314 

0.0304 

^R^vmax^ 

0.0687 

0.1657 

0.0635 

L(m) 

0.1 

0.1 

0.02 

o00(kg/m3) 

0.882 

0.221 

0.882 

0  i 
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The  optical  path  length  (flow  separation  thickness)  was  determined  from  Schlieren 
photographs.  The  density  was  found  by  using  Equation  (2-3).  The  pressure  used 
is  taken  from  the  AKDC  model  atmosphere2^  for  the  given  altitude,  and  the 
temperature  was  taken  as  297°K. 

2 

3.  Results.  Equation  (4-20)  shows  that  a  is  a  function  of  L,  a,  and  the  - 
variance  of  refractive  index  y"^.  Unfortunately  we  do  not  have  a  good  model  for  y2  in 
a  separated  flow.  It  is,  however,  reasonable  to  assume  that  y"2”  increases  with 
freestream  density.  We  now  have  the  problem  of  writing  a2  and  a  as  functions  of 
Poo  and  L  such  that  Equation  (4-61)  is  satisfied  for  each  of  the  three  configurations 
in  Table  II.  The  turbulent  boundary  layer  model  in  which  o2  is  proportional  to 
L2p2  does  not  work.  A  model  which  does  work,  however,  is 


a  *  0.067  L 

2  -12 
a  -  1.2  x  10  A  p  L 


(4-76) 

(4-77) 


All  units  are  MRS.  These  two  equations  along  with  Equation  (4-17)  consititute  our 
model  for  the  average  separated-flow  OTF.  Figure  38  plots  the  resulting  average 
OTF's  (including  both  separated  flow  and  diffraction)  for  the  three  configurations 
and  the  TM.  The  resolution  data  points  Mmy  from  Table  II  are  also  shown. 

As  with  the  MTF  model  for  a  TBL,  there  is  doubt  as  to  the  accuracy  of  our 
model  for  the  average  separated-flow  OTF.  But  again,  it  is  this  writer's  belief 
;  that  the  model  presented  here  is  the  best  that  can  be  found.  The  best  hope  of 
(  achieving  accurate  models  for  the  TBL  and  separated  flow  is  to  take  direct  random- 
flow  MTF  measurements  with  a  shearing  interferometer  under  controlled  conditions. 
The  feasibility  of  such  measurements  has  already  been  demonstrated  by  Kelsall. 

The  Shearing  Interferometer 

The  purpose  of  this  subsection  is  to  present  a  description  of  the  shearing 
interferometer  and  its  use  by  Kelsall  in  aircraft  flight  tests.  Although  several 
complications  make  it  impossible  to  properly  analyze  the  resultant  data,  the 
experiment  is  Important  because  Kelsall  has  demonstrated  the  feasibility  of 
making  direct  measurements  of  the  boundary-flow  MTF. 

Corner-Cube  Shearing  Interferometer .  The  shearing  interferometer  is  an 
interferometer  which  can  be  used  to  directly  measure  MTF's.  Figure  39  is  a 
simplified  diagram  of  a  corner- cube  shearing  interferometer  introduced  by 
Hopkins. 2^  After  passing  through  an  aperture  the  incident  beam  is  divided  by  a 
beam  splitter.  The  waves  are  then  each  reflected  by  a  corner-cube  retroref lector 
and  superimposed  at  the  beam  splitter  with  a  lateral  shear  x.  The  power  produced 
by  the  superimposed  beams  is  recorded.  To  guarantee  that  the  two  wavefronts  are 
in  phase,  one  of  the  corner-cubes  may  be  driven  forward  and  back  over  several 
wavelengths.  The  power  output  will  then  consist  of  an  ac  signal  corresponding 
to  the  overlapping  portion  of  the  wavefronts  and  a  dc  signal  corresponding  to 
the  power  that  would  be  received  if  there  were  no  overlap.  When  the  ac  signal 
is  at  its  maximum,  the  wavefronts  are  exactly  in  phase.  Kelsall  has  shown  that 
the  amplitude  of  the  ac  signal  is  proportional  to  the  MTF,  -r(x/Af),  associated 
with  the  incoming  wavefront.22  The  radiant  output  power  pow  of  the  interferometer 
is  given  by 
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BASIC  DESIGN  OF  A  CORNER-CUBE  SHEARING  INTERFEROMETER 
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pow  =  a  +  8  | t  O^y)  I  cos  [<Kx)  +  upt]  (4-78) 

where  a  and  8  are  constants  (8  ±  a),  $  is  a  phase  term  which  includes  the  OTF 
phase,  and  wp  is  the  angular  frequency  of  the  relative  phase  of  the  two  wavefronts. 
Figure  40  shows  plots  of  power  vs  time  for  various  shear  values,  x/D,  where  D  is 
aperture  diameter.  In  Figure  40(a)  the  wavefront  coming  into  the  interferometer 
is  constant  in  time,  and  in  Figure  40(b)  the  incoming  wavefront  varies  in  time  as 
would  a  wavefront  passing  through  a  TBL  or  separated  flow  region. 

The  average  MTF  for  a  particular  spatial  frequency  is  found  by  taking  the 
average  ac  power  amplitude  for  the  corresponding  value  of  x  and  normalizing  it 
such  that  the  normalized  ac  power  amplitude  at  x  =  0  is  unity. 

Kelsall's  derivation  examines  the  effects  of  finite  bandwidth,  inequality 
of  power  of  the  two  beams,  and  partial  incoherence  caused  by  unequal  beam  paths. 

All  of  these  effects  increase  the  dc  power  level  and  decrease  the  ac  power 
amplitude,  but  the  MTF  is  still  proportional  to  the  ac  power  amplitude.  A 
simpler  derivation  which  does  not  examine  these  effects  but  deals  only  with  the 
basic  principle  of  the  shearing  interferometer  is  presented  in  Appendix  G. 

Kelsall's  Experiment.  The  corner-cube  shearing  interferometer  used  by 
Kelsall  is  described  in  Reference  23.  The  corner-cube  which  moves  forward  and 
back  to  produce  the  phase  shift  between  the  wavefronts  is  driven  by  a  piezo¬ 
electric  device.  Two  identical  glass  plates  are  placed  in  the  return  paths  from 
the  corner-cubes  to  the  beam  splitter.  By  turning  these  plates  simultaneously 
through  a  drum  belt  coupling,  the  shear  x  can  be  adjusted.  These  shear  plates 
and  the  piezo-electric  driver  are  outlined  with  dotted  lines  in  Figure  39. 

Figure  41  shows  the  experimental  arrangement  used  by  Kelsall  to  measure 
the  turbulent  boundary  layer  MTF  of  a  KC-135  aircraft.^  The  source  beam  from  a 
He-Ne  laser  is  expanded  to  50  mm  then  reflected  through  the  aircraft  window  by  a 
plane  mirror  and  a  beam  splitter.  The  beam  is  reflected  back  through  the  window 
by  a  31  mm  mirror  mounted  on  an  external  airfoil  25  cm  from  the  window.  In  this 
way  the  beam  passes  through  the  same  boundary  flow  region  twice.  The  return  beam 
is  then  reduced  by  a  24  mm  aperture  telescope  and  passes  to  the  shearing 
interferometer.  A  removable  air  path  shield  allows  MTF  measurements  to  be  made 
without  the  boundary  flow.  Figure  42  shows  the  MTF's  for  the  double  pass.  The 
measurements  were  taken  at  40,000  ft  altitude  at  normal  cruising  speed. 

We  will  not  attempt  to  analyze  the  results  of  this  experiment,  but  some 
comments  will  be  made.  Obviously  the  MTF  detected  by  the  interferometer  is  not 
the  MTF  of  the  boundary  flow  since  the  beam  makes  two  passes  through  the  flow. 

To  find  the  true  value  of  tr  from  the  double  pass  value  Tr2,  Kelsall  assumed 
that  -ln(t£)  is  proportional  to  the  total  optical  path  length.  Using  our  model 
for  Tr,  Equation  (4-17),  which  is  also  used  by  Kelsall,  the  assumption  implies 
that  o«-  is  proportional  to  the  total  optical  path  length.  The  relationship 
between  Tr  and  Tr  would  be 


T 
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FIGURE  42 

AVERAGE  MTF'S  FOR  DOUBLE  PASS  IN  KELSALL'S  EXPERIMENT6 
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Kelsall's  assumption  would  be  correct  if  the  optical  length  of  the  return  path 
of  the  beam  were  independent  of  the  optical  length  of  the  first  path.  Then  cr? 
for  the  double  pass  would  be  the  variance  of  the  sum  of  two  independent  values 
of  A,  a  Gaussian  variable.  Using  a  law  from  probability  theory,  the  a 2  for  the 
double  pass  would  be  the  sum  of  the  c2  for  the  two  separate  passes.  In  other 
words,  Kelsall's  assumption  would  be  true.  However,  the  optical  lengths  are 
not  independent.  In  fact,  if  the  instrumentation  is  properly  aligned,  the  two 
optical  lengths  are  the  same,  and  the  value  of  A  is  twice  as  large  at  every 
point  for  the  double  pass.  According  to  Equation  (4-8),  ct2  for  the  double  pass 
would  be_four  times  greater  than  foT  a  single  pass.  Therefore  the  relationship 
between  tr  and  should  be 
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R2 


(4-80) 


By  using  Kelsall's  data  and  Equation  (4-80)  one  should  be  able  to  find  tr  for 
the  KC-135  boundary  flow.  Kelsall  also  ran  a  wind  tunnel  test  which  was  to  duplicate 
the  flight  test  conditions.  *  The  results  of  the  wind  tunnel  test  showed  no 
significant  optical  degradation.  Several  factors  may  be  responsible  for  the 
discrepancy  between  the  results  of  the  flight  and  wind  tunnel  tests. 

1.  The  optical  window  was  mounted  flush  to  the  wall  in  the  wind  tunnel 
test;  whereas,  in  the  flight  test  the  window  was  recessed  about  1  inch. 

2.  In  the  flight  test  a  temperature  gradient  of  about  30°F  existed  between 
the  skin  and  the  outside  air. 

3.  It  is  possible  that  during  the  flight  test  warm  air  from  inside  the 
aircraft  escaped  from  openings  in  the  skin  forward  of  the  test  window 
and  mixed  with  the  boundary  flow  air. 

All  three  of  these  factors  would  tend  to  increase  degradation  in  the  flight  test. 

They  also  make  a  proper  analysis  of  the  results  with  regard  to  a  boundary  layer 
MTF  very  difficult  if  not  impossible. 

Although  the  discrepancy  between  the  flight  and  wind  tunnel  test  results 
shows  that  the  flight  test  results  cannot  be  applied  to  a  general  situation,  the 
test  is  important  because  it  shows  the  feasibility  of  making  direct  measurements 
of  boundary  flow  MTF's. 

Until  now  our  discussion  has  been  confined  to  the  time  average  MTF.  New 
fast  shearing  interferometers  now  make  it  possible  to  take  MTF  measurements 
almost  instantaneously  (~  one  millisecond). 

Fast  Shearing  Interferometers.  There  are  two  types  of  fast  shearing  inter¬ 
ferometer;  the  first  was  developed  by  Kelsall25  and  the  second  was  developed  by 
workers  at  MCAIR.  Basic  designs  of  the  two  instruments  are  shown  in  Figure  43. 

The  two  designs  are  similar,  but  because  the  MCAIR  version  is  easier  to  describe, 
it  will  be  discussed  first.  % 
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(a)  MCAIR  Design 
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FIGURE  43 

BASIC  DESIGNS  OF  THE  MCAIR  AND  KELSALL  FAST  SHEARING  INTERFEROMETERS 
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As  shown  in  Che  diagram  of  the  MCAIR  fast  shearing  interferometer,  the 
incident  beam  passes  through  an  aperture  and  is  divided  by  a  beam  splitter.  The 
two  beams  traverse,  in  opposite  directions,  the  triangular  path  formed  by  the 
splitter  and  two  plane  mirrors.  The  beams  are  recombined  at  the  splitter  with 
some  lateral  shear  x.  For  small  angular  changes,  the  shear  is  a  linear  function 
of  the  angle  of  the  beam  splitter.  In  addition,  the  optical  path  difference 
changes  rapidly  with  the  splitter  angle.  These  facts  make  both  the  shear  plates 
and  the  piezoelectric  driver  used  in  the  corner-cube  version  unnecessary.  The 
splitter  is  mounted  through  flexural  pivots  to  the  interferometer  frame  and 
oscillates  from  power  supplied  by  an  electrodynamic  driver.  Fortunately,  an 
oscillation  amplitude  of  only  a  few  degrees  is  necessary  when  the  aperture  diameter 
D  is  one  the  order  of  one  inch.  A  direct  MTF  reading  is  obtained  from  an 
oscilloscope  display  by  using  the  radiant  power  of  the  sheared  beams  and  the 
angular  displacement  of  the  splitter  as  the  vertical  and  horizontal  scope  inputs 
respectively.  Figure  44  shows  two  traces  for  a  rectangular  aperture,  one  with  no 
degradation  and  the  other  with  random  degradation.  Notice  that  each  trace  actually 
produces  four  MTF's,  one  in  each  quadrant. 

In  principle,  Kelsall's  interferometer  shown  in  Figure  43  differs  from  the 
MCAIR  version  only  slightly.  Shear  change  is  produced  by  rotating  a  shear  plate 
in  the  path  between  the  two  mirrors  at  constant  angular  rate  instead  of  oscillating 
the  beam  splitter.  Both  shear  and  optical  path  difference  are  linear  functions  of 
the  shear  plate  angle,  therefore,  the  time  trace  of  the  output  power  yields  the  MTF. 

These  two  fast  shearing  interferometers  provide  a  method  for  measuring 
directly  boundary  flow  MTF's.  The  differences  between  the  two  approaches  are 
subtile.  Kelsall's  instrument  provides  a  constant  shear  rate  in  contrast  to  the 
variable  shear  rate  inherent  in  the  MCAIR  design.  The  MCAIR  approach  should  be 
capable  of  producing  a  faster  MTF  because  of  the  lower  mechanical  motion  required. 
The  faster  the  MTF  is  produced,  the  less  sensitive  it  is  to  vibration  and  acoustic 
noise.  More  Information  on  Kelsall's  Instrument  is  available  in  Reference  25. 

A  detailed  description  of  the  MCAIR  version  will  be  published  in  a  future  MDC 
report. 
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(a)  No  Degradation 


(b)  Random  Degradation 
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FIGURE  44 

SCOPE  TRACE  OF  FAST-SHEARING-INTERFEROMETER  OUTPUT  POWER 
FOR  A  RECTANGULAR  APERTURE 
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V.  Parametric  Analysis 

The  previous  sections  presented  the  optical  quality  of  turbulent  boundary 
layers,  separated  flow,  and  main  flow  in  terms  of  either  the  OTF  or  phase  error  A. 
In  this  section,  the  OTF  s  and  phase  errors  will  be  used  to  compute  the  optical 
effect  of  boundary  flow  on  propagating  and  imaging  systems.  Beam  cross-sections 
and  receiving  apertures  are  circles  of  diameter  D.  The  effect  is  presented  in 
terms  of  normalized  central  power  density  (NCPD)  vs  beam  diameter  for  propagating 
systems  and  in  terms  of  resolution  vs  aperture  diameter  for  Imaging  systems.  The 
following  parameter  values  are  considered. 


Mach  number 

Altitude 

Wavelengths 


0.9 

8000  ft,  40,000  ft 
0.55p,  10. 6\i 


Standard  atmospheric  conditions  at  the  two  altitudes  are  as  follows. 


8000  ft 


40,000  ft 


Density  (kg/ni  )  0.940 

Pressure  (nt/m^)  75,000 

Temperature  (#K)  278 

Speed  of  sound  (m/sec)  334 


0.303 

18,800 

217 

295 


In  addition  to  main  flow,  turbulent  boundary  layers,  and  separated  flow, 
wake  flow  is  also  considered.  It  is  treated  as  a  special  case  of  separated  flow. 
Shock  waves  are  not  considered  because  the  only  quantitative  work  done  here  in 
that  area  con .  irns  beam  direction  only,  not  beam  quality. 

Often  it  is  desirable  to  predict  the  optical  effects  corresponding  to  a 
full  scale  aircraft  by  using  a  scale  model.  The  last  subsection  derives  scaling 
rules  for  achieving  this  objective. 

Propagating  Systems 

The  purpose  of  this  subsection  is  to  present  the  NCPD  as  a  function  of  beam 
diameter  for  main  flow,  TBL's,  separated  flow,  and  wake  flow. 

In  Section  II  it  was  shown  that  the  NCPD  could  be  written  as  a  function  of 
either  the  OTF, 

D 

NCPD  -  f  T(ff)  r  dr,  (5- 


or  the  optical  phase, 


■JW 


.iki(r)r 


Main  Flow.  The  first  step  in  finding  the  NCPD  for  the  main  flow  of  a  sphere 
is  to  find  A.  Equation  (3-15)  was  evaluated  for  ■  0,  that  is,  for  the  beam 
pointed  directly  into  the  flow.  The  sphere  diameter  used  was  one  meter  and  the 
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i  upper  limit  of  the  integral  was  set  at  l  m  50,  which  effectively  is  the  same  as  I  «  •. 
From  the  symmetry  of  the  geometry,  A  is  symmetrical  with  respect  to  the  origin, 

A(x,y)  ■  A(r).  Figure  45  is  a  plot  of  A(r)  for  altitudes  of  8000  ft  and  40,000  ft. 

An  approximation  to  the  main  flow  NCPD  can  be  easily  obtained  by  approximating  A(r) 
with  a  parabola  of  the  form 

AD(r)  -  w  r2.  (5-3) 

Equation  (5-3)  is  plotted  as  a  dotted  line  in  Figure  45  with  A  ■  Aq  at  r  ■  0-2 
meters. 

Inspection  of  Equation  (3-15)  shows  that  for  constant  values  of  ?  *  x/E  and 
n  “  y/R»  A  is  proportional  to  p<„R,  where  p»  is  the  freestream  density^  Equation 
(5-3)  shows  that  for  constant  values  of  r/R,  An  is  proportional  to  w  R  .  There¬ 
fore  if  A(r)  is  approximated  by  AD(r),  then  wR*  is  proportional  to  p„R,  or  w  is 
proportional  to  P<»/R.  If  A(r)  is  set  equal  to  Ac(r)  at  r  ■  0.2  meters,  the  equation 
for  w  is 

w  -  -0.823  x  10~4  Poo/R*  (5-4) 

All  units  are  MRS;  w  has  units  of  meters-1. 

When  a  wavefront  has  a  phase  given  by  Equation  (5-3) ,  it  is  said  to  suffer 
a  defect  of  focus.  The  resulting  NCPD  is  found  by  substituting  Equation  (5-3) 
into  Equation  (5-2) .  The  result  is 

'  HCPD .  m  (5-5) 

(WT/A) 

Equation  (5-5)  is  plotted  in  Figure  46.  The  NCPD  for  main  flow  is  found  by 
substituting  Equation  (5-4)  into  Equation  (5-5).  Figure  47  plots  the  NCPD  as  a 
function  of  the  beam  size  D  for  altitudes  of  8000  ft  and  40,000  ft  and  for 
wavelengths  of  0.55U  and  10. 6p.  The  plots  show  that  the  main~flow  effect  is  not 
serious  for  10. 6p  at  40,000  ft.  For  the  other  conditions,  however,  the  central 
power  may  vanish  completely,  depending  on  the  beam  diameter. 

t 

Turbulent  Boundary  Layer.  The  average  0TF  for  turbulent  flow  is  given  by 
Equations  (4-9)  and  (4-10). 

Uff)  -  T0(jj)  expj-k2c2[l-e"(r/a)2J  J  (5-6) 

With  the  diffraction  0TF  for  a  circular  aperture.  Equation  (2-35),  Equation  (5-6) 
becomes 

T  (fj)  -  |  Icos-1  (§>-(§) /l-(§)2]  expj-k2o2U-e“(t/a)2]},  0  <  r  <  D.  (5-7) 
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FIGURE  45 

PHASE  ERROR  A(r)  FOR  SPHERE  WITH  FORWARD-LOOKING  BEAM 
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For  the  TBL  on  a  flat  plate,  a  and  a  are  found  using  Equations  (4-53)  and  (4-54) 
respectively.  Figure  48  plots  as  a  function  of  D  the  NCPD  found  by  substituting 
Equations  (5-7),  (4-53)  and  (4-54)  into  Equation  (5-1).  The  figure  shows  NCPD 

plots  for  altitudes  of  8000  ft  and  40,000  ft  and  for  wavelengths  of  0.55p  and  10. 6y. 

Distance  from  the  leading  edge  of  the  flat  plate  is  two  meters.  The  plots  show  that 

the  TBL  effect  is  generally  weak.  In  fact  it  is  completely  negligible  at  both 

altitudes  for  the  10. 6p  wavelength.  At  the  shorter  wavelength,  the  NCPD  drops  to 
minimums  of  0.85  and  0.98  for  altitudes  8000  ft  and  40,000  ft  respectively.  The 
fact  that  the  NCPD  remains  constant  at  its  minimum  value  for  all  but  the  smallest 
aperture  sizes  is  a  result  of  the  small  spatial  scale  of  the  wavefront  disturbance. 

Separated  Flow.  The  NCPD  for  separated  flow  is  found  by  again  substituting 
Equation  (5-7)  into  Equation  (5-1).  The  values  of  o2  and  a  are  found  by  using 
Equations  (4-76)  and  (4-77)  respectively.  Figure  49  plots  NCPD  as  a  function  of 
D  for  altitudes  of  8000  ft  and  40,000  ft  and  for  wavelengths  of  0.55P  and  10. 6p. 

The  separation  thickness  is  10  cm.  The  plots  show  the  strong  dependence  of  the 
separated  flow  effect  on  wavelength.  At  10. 6u  the  NCPD  drops  to  only  0.96  and  0.99 
for  altitude  of  8000  ft  and  40,000  ft  respectively.  At  0.55p,  however,  the  NCPD 
drops  to  less  than  0.01  for  both  altitudes.  As  with  the  turbulent-boundary- layer 
NCPD,  the  small  spatial  scale  of  the  wavefront  disturbance  causes  the  separated- 
flow  NCPD  to  be  constant  for  all  but  the  smallest  values  of  D. 


Wake  Flow.  In  this  subsection  the  NCPD  for  the  wake  flow  behind  a  sphere 
one  meter  in  diameter  is  considered.  The  problem  is  treated  as  a  special  case  of 
separated  flow.  Figure  50  shows  the  problem  geometry.  The  angle  between  the 
polar  axis,  which  points  into  the  flow,  and  the  direction  of  propagation  is  Bq. 

For  purposes  of  calculating  the  correlation  length,  a,  the  separation  thickness 
is  taken  to  be  equal  to  the  sphere  radius  R.  Equation  (4-77)  becomes 

a  -  0.067R.  (5-8) 


The  average  optical  path  length  through  the  wake  flow  is  given  by 


L 


1). 


(5-9) 


2 

Since  Equation  (4-20)  implies  that  a  is  proportional  to  the  optical  path  length. 
Equation  (4-76)  becomes 


1.2  x  10 


-12 


P»R(- 


sin  0Q 


-  1). 


(5-10) 


Figure  51  is  a  plot  of  the  resulting  NCPD  as  a  function  of  beam  width  for  altitudes 
of  8000  ft  and  40,000  ft  and  for  wavelengths  of  .55w  and  10. 6p.  The  sphere  radius 
R  is  one-half  meter,  and  the  polar  angle  6o  is  135®.  The  NCPD  for  wake  flow  is 
similar  to  the  NCPD  for  separated  flow.  The  a2  value  is  larger  for  the  wake  flow 
because  of  the  longer  optical  path.  This  larger  value  of  o?  causes  the  NCPD  to 
drop  to  lower  values  for  the  10. 6u  wavelength,  0.85  at  8000  ft  and  0.95  at  40,000 
ft.  On  the  other  hand,  the  value  of  the  correlation  le  ~.h  is  larger  for  the 
wake  flow  than  for  the  separation  flow.  The  result  is  that  the  NCPD  reaches  its 
*1  minimum  value  at  larger  values  of  D. 


FIGURE  SO 

WAKE-FLOW  PROBLEM  GEOMETRY 
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Propagation  Angle  80  ■  135° 
Sphere  Diameter  =  1  meter 
Mach  Number  =  0.9 


FIGURE  51 

NORMALIZED  CENTRAL  POWER  DENSITY  (NCPD)  n  BEAM  DIAMETER  FOR 
THE  WAKE  FLOW  OF  A  SPHERE 
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Suninary.  In  general  boundary  flow  has  only  a  slight  effect  on  10. 6u  radiation.  I 
The  only  exception  Is  the  effect  of  main  flow  at  low  altitudes  in  which  the  NCPD 
may  drop  below  one  percent.  At  a  wavelength  of  0.55p  main  flow,  separated  flow, 
and  wake  flow  can  all  cause  the  NCPD  to  drop  below  one  percent,  depending  on  the 
beam  diameter.  Only  the  TBL  effect  is  not  drastic. 


Imaging  Systems 


A  common  criterion  for  measuring  the  performance  of  an  imaging  system  is 
system  resolution  Vi,aT.  The  units  of  are  inverse  distance,  for  example. 

The  value  of  Vimv  is  found  by  using  Equation  (4-57), 


mm~l  or  lines/mm. 


-  ColT('W)  I  •  0  i  c0  i  (5-ID 

where  |  t (Vmn-r)  I  is  the  MTF  of  the  boundary  flow,  aperture,  and  lens,  IM(v)  is  the 
threshold  modulation  of  the  film  and  human  recorder,  and  Cq  is  the  target  modulation. 

In  an  Idealized  system  in  which  TM(v)  vanishes  and  t(v)  is  diffraction  limited, 
the  resolution  vMT  is  equal  to  the  spatial  frequency  vQ  at  which  the  diffraction 
MTF,  |tq(v) |  vanishes. 


lT0(vmax)l  "  lT0(Vl  “  0  (5_12) 

From  Equation  (2-35)  is  given  by 

r  -  AfvQ  -  D  (5-13)  ) 

or 


where  D  is  aperture  diameter,  A  is  wavelength,  and  f  is  focal  length. 

To  find  VmflY  for  a  real  system,  Co,  t(v),  and  TM(v)  must  all  be  specified. 
However,  there  is  a  method  of  approximating  the  effect  of  boundary  flow  on  image 
resolution  by  using  only  the  MTF  for  the  aperture  and  the  boundary  flow.  This 
method  allows  us  to  examine  the  effect  of  boundary  flow  on  imaging  systems  without 
specifying  a  particular  system.  Let  ?l(v)  and  \>l  be  OTF  and  resolution  respectively 
for  the  system  excluding  boundary  flow.  Let  r(v)  and  vMT  be  the  OTF  and  resolution 
respectively  for  the  system  including  boundary  flow.  Figure  52  shows  plots  of 
Co|tL(v) |  and  CgjT(v)|.  We  assume  that  is  small  enough  so  that  the  IM  is  fairly 
constant  for  v<vL.  This  assumption  is  usually  accurate  for  aerial  reconnaissance 
situations  because  Cg  is  small,  say  0.1  to  0.3.  Two  new  linear  OTF's,  t£'(v)  and 
t~(v),  are  constructed  such  that  the  NCPD's  corresponding  to  tl  and  t£-  are  equal 
and  the  NCPD's  corresponding  to  t  and  are  equal.  The  resolutions  corresponding 
to  and  t"  are  v£'  and  respectively.  The  equations  for  and  x"  are 
respectively 


.£'(v)  -  1  -  -  1  -  ^55-.  0<  .-feih,. 


(5-is;q 


Modulation 


4' 
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FIGURE  S2 

APPROXIMATION  OF  SYSTEM  RESOLUTION  USING  A  LINEAR  OTF 
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I 


and 

x"(v)  -  1  -  v/h-  1-fjj,  0<  v  h 

where  and  h  are  defined  by 

T'£(hL)  -  t"(h)  -  0. 


(5-16) 


(5-17) 


The  values  of  the  NCPD's  are  found  in  terms  of  h^  and  h  by  substituting  Equations 
(5-15)  and  (5-16)  respectively  into  Equation  (5-1). 


XfhL 


NCPDT 


l-4  /  « -  xfe?  ** ■  ! <T>: 


Xfh 


■»  •  $  /"  «  -  fe>  r-r  -  f  <f>)2 


(5-18) 


(5-19) 


NCPD^  is  the  NCPD  corresponding  to  and  *£'. 

The  ratio  vma-r/ vt.  is  approximated  by 
52  we  have 


h 

V 


From  the  geometry  in  Figure 

(5-20) 


Solving  Equations  (5-18)  and  (5-19)  for  hj,  and  h  respectively  and  substituting 
into  Equation  (5-20)  we  arrive  at  the  approximation  formula. 


v 

~L 


v 

j 

VL 


NCPD 

NCPDj/ 


(5-21) 


Let  us  now  work  a  sample  problem  with  the  following  input  parameters. 

Altitude  8000  ft 


Wavelength  .  55 

Focal  Length  152.2  nm 

Aperture  size  9.5  mm  (f/16) 

Target  contrast  2:1 

Type  of 

degradation  10  cm  flow  separation 


Mach  No.  0.9 

Film  Kodak  3414 

Lens  Degradation  None 


First  the  resolution  va „„  is  found  using  Equation  (5-11).  The  TM(v)  for  Kodak  3414 
film  is  given  by  Equation  (4-65).  From  the  target  contrast  of  2:1  and  from  Equation 
(2-50),  Cq  must  satisfy 


*  *  cc 

1  -  C* 


or 


co-  J- 
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(5-22) 


(5-23) 


0  ! 
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The  MTF  is  given  by  Equations  (5-7),  (4-76)  and  (4-77).  For  an  altitude  of  8000 
ft,  p.  is  0.940  nt/m^,  The  value  of  L  is  0.01  m.  Figure  53  plots  both  TM(v) 
and  CqT(v)  as  a  function  of  v.  The  resolution  is  found  to  be  Unnv  »  30  lines/mm. 

The  function  CoTq(v)  is  also  plotted  in  Figure  53.  The  corresponding  diffraction 
limit  resolution  is  vp  ■  84  lines/am.  To  use  the  approximation  formula.  Equation 
(5-21)  we  set  vj,  *vp  -  84  lines/mm  and  NCPDl  "  1  since  there  is  no  lens  degradation. 
From  Figure  49  the  value  of  NCPD  is  0.115.  Using  the  approximation 


JWLm  »»  ^  <5"24> 

the  system  resolution  MimT  becomes  28  lines/mm.  We  see  that  for  this  case 
our  approximation  for  ^  is  good. 

Often  it  is  convenient  to  express  the  performance  of  an  imaging  system  in 
terms  of  angular  resolution  AR  or  ground  resolution  GR.  Figure  54  shows  geometrically 
how  AR  and  GR  are  related  to  Nfaax-  The  relationship  may  be  expressed  as  follows. 


AR 


1 

^vmax 


(5-25) 


GR 


SR-AR 


SR 


fv, 


max 


(5-26) 


The  term  SR  is  the  slant  range  from  the  ground  target  to  the  imaging  system.  In 
the  idealized  imaging  system  corresponding  to  Equation  (5-14),  AR  and  GR  become 
respectively 


AR  -  ^  (5-27) 

and 

GR  -  SR  p  (5-28) 

Figures  55  and  56  show  vMT  and  the  corresponding  values  of  AR  and  GR  as  functions 
of  D  for  the  different  types  of  boundary  flow.  In  Figure  55,  altitude  is  8000  ft, 
in  Figure  56,  40,000  ft.  Also  shown  are  the  diffraction-limit  resolutions  vj). 

The  system  resolution  Vim,  is  calculated  using  the  approximation  formula, 

Equation  (5-24).  Values  of  wavelength,  focal  length,  target  contrast,  mach  number, 
and  lens  degradation  are  the  same  as  those  in  the  sample  problem.  Conditions  for 
each  type  of  flow  are  the  same  as  those  used  in  the  Propagation  Systems  subsection. 
In  calculating  GR,  vertical  imaging  was  assumed,  so  that  slant  range  equals  the 
altitude. 


issri 
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FIGURE  54 

RELATIONSHIP  BETWEEN  RESOLUTION  »MAX,  ANGULAR  RESOLUTION  AR, 
AND  GROUND  RESOLUTION  GR 
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Scaling 

This  subseccion  derives  Che  scaling  rules  for  predicting  Che  boundary-flow 
optical  effects  of  a  body  from  a  similar  body  of  different  scale.  The  derivation 
uses  the  aerodynamic  scaling  rules  from  Section  II. 

We  specify  here  that  when  Che  body  is  scaled ,  the  aperture  area  or  beam  cross- 
section  area  I  is  scaled  by  the  same  factor.  For  the  moment  we  assume  that  all 
linear  dimensions  of  the  boundary  flow  are  also  scaled  by  the  same  factor.  Let  the 
aerodynamic  system  containing  the  original  body  be  defined  as  System  1,  and  let 
the  system  containing  the  scaled  body  be  System  2.  Let  S  be  the  scale  of  Body  2 
relative  to  Body  1.  If  the  two  bodies  are  the  same  size,  S  is  one.  Given  the 
temperature  Ti,  Mach  number  Mi,  density  pj_,  and  wavelength  X^  of  System  1  and  the 
scale  S,  we  wish  to  find  the  appropriate  temperature  T£,  Mach  number  M£,  density 
P2>  and  wavelength  X2  for  System  2  such  that  the  two  systems  produce  similar 
optical  effects.  Similar  optical  effects  are  produced  if  the  pupil  functions  are 
similar. 


G^x.y)  -  G2(Sx,Sy)  (5-29) 

Since  the  amplitude  part  of  the  pupil  function  is  constant,  Equation  (5-29)  is 
satisfied  if  the  phases  are  similar, 

k1A1(x,y)  -  k2A2(Sx,Sy)  (5-30) 

or 

A1(x,y)/X1  -  A2(Sx,Sy)/X2.  (5-31) 

Given  Tj_,  Mj_,  pj_,  Xj_,  and  S  the  problem  is  to  find  T2,  M2,  P2*  and  X2  such  that 
Equation  (5-31)  is  satisfied.  The  equation  relating  the  two  areas,  £1  and  £2* 

£2  -  S2!  (5-32) 

and  Equation  (2-28)  show  that  if  Equation  (5-31)  is  satisfied,  then  the  OTF's  of 
the  two  systems  are  similar. 


(5-33) 


Likewise  Equation  (2-51)  show  that  if  Equation  (5-31)  is  satisfied,  then  the 
NCPD's  of  the  two  systems  are  equal. 


NCPD^  -  NCPD2 


(5-34) 
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General  Case.  Rules  are  now  derived  for  finding  appropriate  values  of  T2,  M2, 
P2»  and  ^2*  These  rules  may  be  applied  to  all  types  of  boundary  flow. 

The  relationship  between  A  and  density  p  is  given  by  Equation  (3-3). 

L 

A(x,y)  ■  K  f  [p(x, y,z)  -  p(0,0,z)]  dz  (5-35) 

0 


Substituting  Equation  (5-35)  into  Equation  (5-31)  yields 


Ll  SL2 

“  J"  tP2(x»y»z)  ■*  p(0,0,z)]  dz  *  J  fp2(Sx,Sy,z)  -  P2(0,0,z)]  dz 


(5-36) 


or 


L  I  w 

1  ^1 
»1  f  pl(x»y*z)  ■  SP.,  f  P2(Sx,Sy,Sz)  -  p  (0,0, Sz) 

—  J  - - - dz  -  — —  J  - - -  dz 


(5-37) 


10  “1  2  0  p®2 
where  p^  and  pm2  are  the  freestream  densities. 

Examination  of  Equation  (5-37)  shows  that  if  the  equations 


Pool  SPm2 

— ».*  m  .  .  Jr 

A1  *2 


(5-38) 


and 


P2(x,y,z)  p2(Sx,Sy,Sz) 


J«1 


(5-39) 


“2 


are  satisfied,  then  Equation  (5-31)  is  satisfied.  Equations  (5-38)  and  (5-39)  are, 
therefore,  the  general  optical  scaling  rules.  Generally  p„2  and  A  2  can  be  adjusted 
so  that  Equation  (5-38)  holds.  Equation  (5-39)  will  hold  if  the  aerodynamic  flows 
for  the  two  systems  are  similar.  Section  XI  shows  that  similarity  of  the  flows 
is  guaranteed  if  the  Mach  numbers  and  Reynolds  numbers  are  equal. 


M 


1 


M 


2 


(5-40) 

(5-41) 
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For  equal  temperatures 


Tx  -  T2  (5-42) 

and  Mach  numbers  ,  Equation  (5-41)  becomes 

Pi-Spj.  (5-43) 

Substituting  Equation  (5-43)  with  freestream  values  of  p-i  and  p?  Into  Equation 
(5-38)  yields 


X 


1 


2* 


(5-44) 


In  terms  of  freestream  values  the  general  optical  scaling  laws  may  be  summarized 
as  follows 


T»1  ”  t®2*  M»l  "  **<■2*  p«i  "  Sp»2  »  *i  *  ^2 


(5-45) 


If  Equations  (5-45)  are  satisfied,  the  two  systems  produce  similar  flows  and 
similar  optical  effects  and  Equations  (5-31),  (5-33),  and  (5-34)  all  hold. 

In  scale  model  testing  It  Is  often  difficult  or  impossible  to  meet  the 
requirement  for  density.  Approximately  similar  optical  effects  can  be  produced, 
however,  even  if  the  requirements  in  Equation  (5-45)  are  relaxed  as  shown  below. 

Main  Flow.  We  have  shown  that  if  Equations  (5-38)  and  (5-39)  are  satisfied 
the  two  systems  produce  similar  optical  effects.  Using  Equation  (3-4),  which  is 
based  on  the  Karman-Tsien  approximation  for  subsonic  main  flow,  we  can  now  show 
that  Equation  (5-39)  Is  automatically  satisfied  for  subsonic  main  flow  if  free¬ 
stream  temperatures  and  Mach  numbers  for  the  two  systems  are  equal. 

Substituting  Equations  (3-5),  (3-6),  and  (3-9)  into  Equation  (3-4)  yields 


2(y-l)M2  +  4  -  V2(x,y,z)/a2 
®  00 

2(y-l)M2  +  4  +  V2(x,y,z)/a* 
00  00 


(5-46) 


where  V  is  the  potential-flow  velocity  field.  It  is  a  property  of  potential  flow 
that  similar  bodies  of  different  scales  produce  similar  velocity  fields 


VA(x,y,z) 


V2(Sx,Sy,Sz) 


(5-47) 
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t.  The  U's  are  freestream  speeds.  If  the  freestream  Mach  numbers  and  temperatures 
are  equal,  then  the  freestream  speeds  are  equal. 


Ul-°2 


(5-48) 


Equations  (5-46),  (5-47),  and  (5-48)  show  that  Equation  (5-39)  is  satisfied.  The 
optical  scaling  rules  for  subsonic  main  flow  may  be  summarized  as  follows 


t-2*  m»i  ■  m-2‘  r 


(5-49) 


Turbulent  Boundary  Layers.  Although  we  have  no  equations  for  A  or  p  in  a 
random  flow,  we  do  have  equations  for  the  average  OTF.  For  this  reason  Equation 
(5-33)  Is  used  as  the  similarity  condition  to  be  met  rather  than  Equation  (5-31) . 

Assuming  that  freestream  temperatures  and  Mach  numbers  are  equal  for  the  two 
systems.  Equations  (2-34),  (4-9) ,  (4-16) ,  (4-53),  and  (4-54)  yield 


2/7..  12/7 


tqi  and  T02  are  the  diffractio 
respectively.  Note  that  the  T q 
1,  the  following  equations  must 
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These  equations  imply 


“  Sp-2 


(5-54) 


(5-55) 


These  are  the  optical  scaling  rules  for  the  general  case.  However,  it  is 
possible  to  relax  Equations  (5-54)  and  (5-55)  and  still  obtain  approximately 
similar  optical  effects  if  the  wavelength  is  properly  adjusted.  Equation  (5-51) 
shows  that  the  exponential  term  varies  according  to  the  two-sevenths  power  of 
p^S  And  according  to  the  second  power  of  x.  Equation  (5-53)  can,  therefore,  be 
ignored  without  too  much  error  provided  that  Sp*^  is  within  a  factor  of  two  or  three 
of  pool*  Then  only  Equation  (5-52)  must  be  satisfied.  This  equation,  however,  is 
approximately  the  same  as  Equation  (5-38).  The  following  scaling  rules  may, 
therefore,  be  used  to  obtain  approximately  similar  optical  effects  for  a  TBL. 


P«1  P“2 

T-1  ‘  T-2‘  M-1  “  “-2*  if  "  S  17 

These  are  the  same  equations  that  were  derived  for  subsonic  main  flow. 


(5-56) 


Separated  Flow.  The  rules  for  seperated  flow  are  found  in  the  same  way  as 
the  rules  for  the  TBL. 

Assuming  that  Mach  numbers  for  both  systems  are  in  the  0.8  -  0.95  range. 
Equations  (2-34),  (4-9),  (4-16),  (4-76),  and  (4-77)  yield 


fi(w^)-  Toi(xiV^r) *xp j-c5  xyf  -  exp(-  xc  +  l )] 


(5-57) 


/_§*_  _sz_\  .  T  /_s*sz_\ 

^2VX2f2j  02^2f2’X2f2; 


exp  -  c. 


/  S2x2  +  S2v2  VI 


(5-58) 
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where  C5  and  %  are  constants.  For  these  two  average  OTF's  to  be  equal  the 
following  equation  must  be  satisfied. 


The  scaling  rules  for  separated  flow  may  be  summarized  as  follows. 


0.8  <  M  ., 

—  o>l 


M 


°°2 


£  0.95, 


(5-59) 


(5-60) 


Since  wake  flow  is  a  special  case  of  separated  flow,  these  rules  also  apply  to  wake 
flow. 


Summary  and  Comnents .  Figure  57  summarizes  the  optical  similarity  rules  for 
the  general,  main-flow,  TBL,  and  separated-flow  cases.  It  must  be  remembered 
that  the  rules  for  main  flow,  TBL,  and  separated  flow  were  derived  from  models 
which  are  only  approximately  correct.  Even  so,  all  of  the  scaling  rules  are 
consistant,  since  satisfaction  of  the  general  rules  implies  satisfaction  of  the 
other  scaling  rules. 

The  rules  for  most  of  the  cases  include  the  equality  of  the  freestream 
temperatures.  It  should  be  noted  that  even  if  these  temperatures  vary  from  each 
other  by  several  degrees,  they  may  differ  by  only  a  few  percent  on  an  absolute 
scale. 
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Appendix  A 

Notations.  Symbols,  and  Important  Equations 


This  appendix  is  divided  into  £ive  parts: 
Acronyms 

Abbreviations  for  Units 
Notation  for  Operations 
Symbols 

Important  Equations 


Acruuyuia 

MTF 

Modulation  Transfer  Function 

NCPD 

Normalized  Central  Power  Density 

OTF 

Optical  Transfer  Function 

TM 

Threshold  Modulation 

TBL 

Turbulent  Boundary  Layer 

Abbreviations  for  Units 

°K 

Degrees  Kelvin 

kg 

Kilogram 

in 

Meter 

mm 

Millimeter 

nt 

Newton 

u 

Micron 

Notation  for  Operations 

E 

Expected  value 

e 

Is  an  element  of 

i 

Is  not  an  element  of 

* 

(superscript)  Complex  conjugate 

— 

(overbar)  Tima  average 

) 
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Symbols 

A 

AR 

B 


C 


N 


NCPD 

HCPDj^ 


Optical-wave  amplitude 
Angular  resolution 

Constant  (Equations  (4-67)  and  (4-68)) 

Radiant  power  emerging  from  an  aperture 
Target  modulation 
Image  modulation 
Fluctuating  pressure  coefficient 
Diameter  of  circular  aperture 
Dimension  of  rectangular  aperture 
Dimension  of  rectangular  aperture 
Optical-wave  function 
Pupil  function 

Diffraction  -  limit  pupil  function 
Pupil  function  for  System  1 
Pupil  function  for  System  2 
Ground  resolution 

Working  variable  (Equations  (3-7)  and  (3-8)) 

Image  irradiance  or  power  density  in  the  focal  plane 
Diffraction  limit  of  I 

Zero-order  Bessel  function  of  the  first  kind 
First-order  Bessel  function  of  the  first  kind 
ith-order  Bessel  function  of  the  first  kind 
Gladstone  -  Dale  constant 
Physical  length  of  optical  path 
Physical  length  of  optical  path  for  System  1 
Mach  number 

Mach  number  for  System  1 
Mach  number  for  System  2 

Free-stream  Mach  number 
Free-stream  Mach  number  for  System  1 
Free-stream  Mach  number  for  System  2 
Refractive  modulus 
Normalized  central  power  density 

Normalized  central  power  density  with  no  boundary-flow  effect 
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Normalized  central  power  denaity  for  System  1 
Normalized  central  power  density  for  System  2 
Object  radiance 
Pressure 

Water-vapor  pressure 
Freestream  pressure 

Normalized  output  power  for  Stine-Winovich  experiment 

Radial  coordinate  (Vx^  +  ) 

Radius  of  sphere 
Reynolds  number 
Scale 

Fourier  transform  of  image  radiance 
Fourier  transform  of  object  lrradiance 
Slant  range 
Temperature 

Temperature  for  System  1 
Temperature  for  System  2 
Freestream  temperature 
Freestream  temperature  for  System  1 
Freestream  temperature  for  System  2 
Threshold  Modulation 

Complex  amplitude  of  optical  wave  in  the  focal  plane 

Freestream  airspeed 

Freestream  airspeed  for  System  1 

Freestream  airspeed  for  System  2 

Speed  of  fluid  (air) 

Air  speed  for  System  1 
Air  speed  for  System  2 
Optical  wavefront 
Characteristic  linear  dimension 
Distance  from  leading  edge  of  a  flat  plate 
Cartesian  coordinate  in  the  focal  plane 
Cartesian  coordinate  in  the  object  plane 
Cartesian  coordinate  in  the  focal  plane 
Cartesian  coordinate  in  the  object  plane 
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a 

a 


c 


f 

fi 

f2 

g 

h 

\ 

k 


n 

nl 

n2 

P 

pow 

q 


Speed  of  sound  in  a  fluid  (air) 

Correlation  length  of  wave  phase  A 
Stagnation  or  reservoir  speed  of  Bound 
Criti-al  speed  of  sound 

Correlation  length  of  wave  phase  A  in  the  x-direction 

Correlation  length  of  wave  phase  A  in  the  y-direction 

Frees tr earn  speed  of  sound 

Constant  in  threshold  modulation  model 

Constant  in  threshold  modulation  model 

Speed  of  light 

Working  variable  (Equations  (3-15)  and  (3-16)) 

Working  variable  (Equations  (3-15)  and  (3-17)) 

Constant  (Equations  (5-50)  and  (5-51)) 

Constant  (Equations  (5-50)  and  (5-51)) 

Constant  (Equations  (5-57)  and  (5-58)) 

Constant  (Equations  (5-57)  and  (5-58)) 

Specific  heat  at  constant  pressure 
Specific  heat  at  constant  volume 
Focal  length 

Focal  length  for  System  1 

Focal  length  for  System  2 

Acceleration  of  gravity 

Spatial  frequency  at  which  x"  vanishes 

Spatial  frequency  at  which  x"  vanishes 

L 

Wave  number  (2ir/X) 

Wave  number  for  System  1 
Wave  number  for  System  2 

Physical  length  of  optical  path  normalized  by  sphere  radius 
Integer 

Refractive  index 

Refractive  index  ahead  of  shock  wave 
Refractive  index  behind  shock  wave 
Point  spread  function 

Radiant  output  power  of  a  shearing  interferometer 
Working  variable  (Equations  (3-15)  and  (3-18)) 
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Radial 


coordinate  0 lx2  +  y2  ) 


Radius  of  source  aperture  is  Stine-Wlnovich  experiment 
Radius  of  collection  aperture  in  Stine-Winovich  experiment 
Distance  along  optical  path 
Time 

Fluid  velocity  component  in  the  x-direction 

Dummy  variable  in  Equation  (4-49) 

Horizontal  velocity  component  ahead  of  shock  wave 

Horizontal  velocity  component  behind  shock  wave 

Normalized  value  of  u^  (u^/ac) 

Normalized  value  of  Uj  (Uj/a^ 

Fluid  velocity  component  in  the  y-direction 

Vertical  velocity  component  behind  shock  wave 

Normalized  value  of  v.  (v./a  ) 

2  2  C 

Spherical  abberation  parameter  (Equation  (5-3)) 

Working  variable  (Equations  (3-7)  and  (3-10)) 

Cartesian  coordinate 
Cartesian  coordinate 
Cartesian  coordinate 
Cartesian  coordinate 
Cartesian  coordinate 
Cartesian  coordinate 
Optical-wave  phase 

Optical-wave  phase  in  terms  of  normalized  coordinates  C  and  n 

Optical-wave  phase  for  System  1 

Optical-wave  phase  for  System  2 

Optical-wave  phase  for  defocused  system 

Variation  of  pressure  from  its  mean  (P-P) 

Refractive  index  change  across  shock  wave  (nj-n^) 

Variation  of  density  from  its  mean  (p-p) 

Density  change  across  a  shock  wave  ^-p^) 

Refraction  angle  for  a  beam  crossing  a  shock  wave 
Aperture  area  or  beam  cross-section  area 
Aperture  area  or  beam  cross-section  area  for  System  1 
Aperture  area  or  beam  cross-section  area  for  System  2 
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«► 

4  * 

a 

B 

B 

Y 

6 

6 

6* 

C 

n 

e 

e 


e 

max 


5 


P 

p0 

C  01 


Covariance  function  of  wave  phase  A 

The  dc  component  of  pow 

The  ac  amplitude  of  pow 

Shock  wave  angle 

Ratio  of  specific  heats  (c^/c^) 

Delta  function 

Boundary  layer  thickness 

Boundary  layer  displacement  thickness 

Cartesian  coordinate  x  normalized  by  sphere  radius 

Cartesian  coordinate  y  normalized  by  sphere  radius 

Polar  coordinate  (Equations  (3-15),  (3-18),  and  (3-19)) 

Ramp  angle  in  plane  oblique  shock  wave  problem 

Maximum  ramp  angle  capable  of  producing  a  plane  oblique  shock  wave 

Angle  between  polar  axis  (direction  pointing  into  flow)  and  direction 
of  propagation 

Coefficient  of  thermal  conductivity 
Optical  wavelength 
Optical  wavelength  in  System  1 
Optical  wavelength  in  System  2 

Variation  of  refractive  index  from  its  mean  (n-YO 
Coefficient  of  viscosity 
Freestream  coefficient  of  viscosity 
Spatial  frequency 

Spatial  frequency  at  which  OTF  vanishes 
Spatial  frequency  in  the  x-direction 
Spatial  frequency  in  the  y-direction 
System  resolution 

System  resolution  with  diffraction  OTF  and  nonzero  TM 
Resolution  corresponding  to 
Resolution  corresponding  to  t" 

Resolution  corresponding  to  t^" 

Cartesian  coordinate  z  normalized  by  sphere  radius 
Density 

Stagnation  or  reservoir  density 
Density  ahead  of  shock  wave 


issr 
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Density  for  System  1 
Density  behind  shock  wave 
Density  for  System  2 
Frees tream  density 
Freestream  density  for  System  1 
Frees tream  density  for  System  2 
Standard  deviation  of  wave  phase  A 
Variance  of  wave  phase  A  (Square  of  a) 

Optical  transfer  function  (OTF) 

Diffraction-limit  OTF 
OTF  for  System  1 
OTF  for  System  2 

Diffraction-limit  OTF  for  System  1 

Diffraction-limit  OTF  for  System  2 

OTF  for  controlled  wavefront  aberrations 

OTF  of  system  with  no  boundary-flow  effect 

OTF  degradation  caused  by  random  wavefront  aberration  (t/t  ) 

V 

Measured  value  of  TR  for  double  pass  through  boundary  layer  in  Kelsall  tee‘. 

Linear  OTF  with  same  NCPD  as  t 

Linear  OTF  with  same  NCPD  as  xT 

L 

Minimum  value  of  T_ 

fv 

Phase  factor 

Angle  between  beam  crossing  a  shock  wave  and  the  normal  to  the  shock 
wave 

Normalized  radial  coordinate  (Vc2  +  n2  +  52  ) 

Angular  frequency  of  optical  wave  (2ircA) 

Angular  frequency  of  relative  phase  of  the  two  wavefronts  in  a 
shearing  interferometer 
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t 


Important  Equations 

All  units  are  MKS  unless  otherwise  noted.  Optical  scaling  rules  are  listed 
in  Figure  57. 

Relationships  Between  the  Aerodynamics  Parameters  P ,  T,  x»  and  a  (Speed  of 

Sound) . 


P  •  287  p  T  (2-3) 

Y  “  1.4 

a  -  /287yT  (2-53) 

Relationships  Between  A_,  n,  K,  and 

n  -  1  +  Kp  (2-5) 

K  -  0.223  X  10-3(1  +  7,52  - )  n3/kg  (2  . 

X- 

K  =  0.223  X  10”3  (2-7) 

L 

A(x,y)  -  J  n(x,y,z)  dz  +  mX  (2-13) 

0 


If  A  is  normalized  according  to  A(0,0)  -  0,  then  the  equation  for  A  is 

L 

A(x,y)  -  f  Tn(x  ,y,z)  -  n(0,0,z>]  dz.  (3-2) 

0 


Likewise, 


L  L 

A(x,y)  • j  [n(x,y, z)  -  n]  dz  ■ J  p(x,y,z)  dz 
0  0 


(E-5) 


Is  used  if  A  is  normalized  according  to  A(x,y)  ■  0 
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then 


NCPD  -  3-U  -  -|os(^  IAHm 
(kwD  /4) 

Wave-Phase  Equations  for  Subsonic  Main  Flow. 

L 

A(x  v)  -  H  f  —  Si S&iXtll  _  L.-.JL  (0, OjJ.)  d_ 
u,yj  V  1  +  w'(x,y,z)  1  +  w'(0,0,z)  dz 


H  -  K  p.  ( 1  +  ^  m3 


1 

Y-l 


2(y-1)M  2  +  4 


The  following  equations  apply  to  a  sphere  of  radius  R. 

C  -  x/R  n  -  y/R  5  -  a/R 

*  -  C2  +  n2  +  £2  £  -  L/R 


-  c,  -  q(0,t|>)  c.  -  q  (6n,e) 
A(x,y)  -  A'(c.n)  -  *.Rc2J  ^  *  c7+-q  (e0,O  de 


Cl  -  [2(y-1)Mw2  +  4]/M_2 


c2  “  <1  +  ^  M.2)^1 


q(e,if0  -  1  +  (1  -  3  cos20)/i p3  +  1/4(1  +  3  cos20)/*6 


cos0  ■  (Ccos0q  -  nsin0Q>/i(/ 


(5-5) 

(3-7) 

(3-8) 

(3-10) 

(3-13) 

(3-14) 

(3-15) 

(3-16) 

(3-17) 

(3-18) 

(3-19) 
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Diffraction  Angle  for  11  Beam  Crossing  a  Plane  Shock  Wave. 


4$  -  -2 KPX 


2  2 

H^sinB  -  1 
(Y-l)M(io2sin26  +  2 


Can  4 


_  2 

Equation  for  _t  ,  a  .  and  a^  for  Random  Plow. 


Tflt*If)  "  T0<Xf  Xf>TRCxf  if) 


2_2r,  _-(r/a)i 


For  a  TBL 


TR(Xf*Xf)  ’  tR  (Xf>  “  “P{-k  0  I1’* 


a2  -  4.4  x  10*1262  _2K  2 


a  -  0.086 


For  a  TBL  on  a  flat  plate 


a2  -  4.0*  ,  10-160.12/7T.1/1V6/V2/7 
.  -  7.66  ,  lO-V^V^V1'7*6'7 


For  a  separated  flow 


2  -12 
a  -  1.2  X  10  PL 


a  -  0.067L 


For  a  wake  flow  behind  a  sphere  of  radius  R 


a2  •  1.7  X  -  1) 


a  -  0.067R 


auatlona  for  Imaging  Svst« 


Vx  "  Xf 


y  ^ 


(3-39) 


(4-16) 


(4-52) 


(4-50) 


(4-53) 


(4-54) 


(4-76) 


(4-77) 


(5-10) 


(5-8) 


(2-25) 


(4-18)  0 
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Threshold  modulation  model 

IM(v)  ■  bQ  +  b1  v2 


For  Kodak  3414  film 

TM(v)  -  0.3  +  3.22  X  10_6\>2 
where  v  is  In  mn  ^  or  lines/mm. 

Resolution  condition 

TM(v  )  •  C Jt(v  )  I 

max  0  max 


Resolution  approximation 


(4-63) 


(4-65) 


(4-60) 


Jncpd 

1NCPDl 


(5-21) 


Relationships  between  system  resolution,  angular  resolution  and  ground  resolution 


AR 


(5-25) 


SR 

GR  -  SR  •  AR  -  f  *  (5-26) 

v 

max 


For  an  idealized  system  in  which  TM(v)  vanishes  and  t(v)  is  diffraction  limited. 


V  ■ 
max 

D 

m  - 

Xf 

(5-14) 

AR  - 

X 

D 

(5-27) 

GR  - 

“b 

(5-28) 
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Appendix  B 

Derivation  of  A(x.y)  for  the  Main-Flow  Sphere  Problem 
Equation#  (3-7)  through  (3-11)  give  the  general  solution  for  A(x,y) . 

A(x,y)  .  »/  ^-;4g;jr;g  -  [  ;  < 

h  -  b„(i  +  m_2)’'-1  ( 

M  -  U/a  < 

00  00 

w'(x,y,z)  -  V*(l 'J'.PltL  ( 

2(y-l)M J  +  4 

V2(x,y,z)  <  a^2  < 

The  purpose  of  this  appendix  is  to  find  V  for  a  sphere  (of  radius  R  from  which 
A(x,y)  may  then  be  found.  The  problem  geometry  is  shown  in  Figure  15.  The 
quantity  V  is  the  air  velocity  calculated  from  potential  flow  theory. 

The  velocity  potential  Y  must  satisfy  the  Laplace  equation, 

V2Y  -  0,  ( 


(B-3) 


(B-4) 


(B-5) 


(B-6) 


and  the  problem  boundary  conditions.  The  velocity  V  is  the  gradient  of  Y.  The 
radial  and  polar  velocity  components,  Vr  and  Vg  respectively,  are  given  in  terms 
of  Y  by 

Vr  ‘  £  «■ 


V  „  A 
0  r  66  ' 


(B-8) 
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From  che  symnetry  of  che  problem  geometry  the  azimuthal  component  V$  is  zero. 

The  boundary  conditions  are  as  follows. 

1.  For  r  »  R 

V  -  -U  cos6  V  -  U  sine  (B-9) 

r  w 

2.  For  r  •  R 


V  -  0 
r 


For  our  problem  (V^  -  0)  Equation  (B-6)  becomes 

1  <5/2  ,  1  6/  J  c  6<r\  n 


(B-ll) 


An  expression  for  V  is  needed  such  that  the  boundary  conditions  and  Equation 
(B-ll)  are  satisfied.  Such  an  expression  is 


cose. 


The  velocity  components  are 


and 


Vr  -  -0(1  -  R3/r3)  cose 


U(1  +  Sl_  j  sine. 
V  2r3/ 


(B-12) 


(B-13) 


(B-14) 


In  order  to  calculate  the  Integral  in  Equation  (B-l),  the  coordinates  r  and  e 
must  be  related  to  the  coordinates  xf  y  and  z  in  Figure  15. 


2  2  2  2 

r  ■  x  +  yz  +  z-  (B-15) 

An  expression  for  e  in  terms  of  x,  y  and  z  is  found  by  taking  a  vector  from  the 
sphere  center  to  the  point  (x,ytz),  end  then  expressing  the  projection  of  this 
vector  on  the  polar  axis  in  both  coordinate  systems.  Equating  the  two  expressions 
yields 


r  cose  ■  z  cos6q  -  y  sin6g. 


(B-16) 
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For  simplicity  ths  following  normalized  coordinates  are  used  throughout  the 
rest  of  this  appendix. 


4  ■  x/R  n  -  y/R  C  -  z/R 

4>2  -  r2/R2  -  ;2  +  n2  +  52  i  •  L/R 
The  normalized  velocity-square  function,  q(0,i|<)  ■  becomes 


(B-17) 


V  ^  +  V  2 

q(0.< l>)  m  — — -  1  +  (1  “  3  cos2e)/i(i3  +  1/4(1  +  3  cos2©)/*6,  (B-18) 


and  Equation  (B-16)  may  be  written 


Equation  (B-4)  becomes 


where 


cos0  -  (£  cos0q  -  ti  sin0Q)/i|i. 


w"(x,y,z)  -  q(0/*)/c. 


c2  -  t2(y-l)Mw2  +  41/H^2. 


(B-19) 


(B-20) 


(B-21) 


Finally  Equation  (B-l)  becomes 


-  c  -  q(0,*)  c  -  q(0  ,5) 

Mx.rt  -  *•<«.,)  -  ft..,  f  q7e,.?  -  « 


(B-22) 


where 


C2  “  (1  +  ^  M-2)Y’1 


(B-23) 


Notice  that  the  lower  limit  is  now  one  instead  of  zero.  It  is  assumed  that  the 
beam  is  plane  parallel  at  z  ■  R  or  5*  1,  and  the  Integral  between  £  ■  0  and 
5  ■  1,  therefore,  vanishes. 


To  satisfy  the  Inequality  (B-5), 

<  /l/q<0,*) 


(B-24) 


) 


V 
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must  be  satisfied  for  all  q(6,i|/),  0  <_  0  <_  it  and  >_  1.  Because  q  is  symmetric 
about  0  ■  it/2,  the  limits  on  0  may  be  reduced  to  0  <  0  <  ir/2.  It  is  necessary  to 
find  the  upper  limit  of  q.  Suppose  that  the  iJT^  term  in  Equation  (B-18)  is  non- 
positive.  The  maximum  value  of  the  term  is  zero,  the  nunHimim  value  of  the 
term  is  one,  and  therefore  q  must  be  less  than  or  equal  to  two.  However, 
q(7r/2,l)  is  equal  to  2  1/4.  The  i)f3  term  must,  therefore,  be  positive  when  q  is 
at  its  maximum.  Since  all  terms  in  Equation  (B-18)  are  positive,  tp  must  equal  its 
minimum  value  of  one  for  q  to  equal  its  maximum.  The  value  of  0  which  maximizes  q 
is  found  by  differentiating  q(0,l)  and  setting  it  equal  to  zero. 


-  ~  COS0  Sin0  -  0  (B-25) 

d0  l 

The  pertinent  solutions  are  0  ■  0  and  0  “  tt/2.  If  6  ■  0,  however,  the  term  is 
negative.  The  maximum  value  of  q(0,4>)  is,  therefore,  given  by 


®ax[q(0,ijj)]  -  q(ir/2,l)  -  9/4. 
The  Inequality  (B-24)  becomes 

<  /l/q"(0  1  ^1/ q  (Tr/ 2 , 1)  - 


(B-26) 

(B-27) 
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Appandix  C 

Derivation  of  the  Density  Change  Across  a,  Plane  Oblique  Shock  Wave 
When  the  Shock-Wave  Angle  8,  is.  Known 

Reference  11  gives  the  following  conservation  equations  for  flow  across  a 
plane  oblique  shock  wave. 

Mass 


Posing  -  P2(u2sinB  -  v2cos8)  (C-l) 

Momentum  normal  to  shock  wave 

P^  +  P^u^sin^B  ■  P2  +  P2(u2sinS  -  v2cosB)2  (C-2) 

Momentum  parallel  to  shock  wave 
2 

Pj^  sinBcosB  -  P2(u2sin8  -  v2cosB) (u2cosB+v2sin8)  (C-3) 

Energy 

P  P 

1/2  "i2  ♦  w  ^  ■  l/2(u22  +  v22>  +  A  (c-4) 

The  symbols  u,  v,  p,  and  P  represent  respectively  the  horizontal  and  vertical 
velocity  components,  density,  and  pressure.  Subscripts  1  and  2  denote  respectively 
values  upstream  and  downstream  of  the  shock  wave.  Equations  (C-l)  through  (C-4) 
are  a  system  of  four  equations  with  four  unknowns,  u2,  v2  >  p2 ’  811(1  °2'  The  problem 
here  is  to  find  p2.  This  task  is  accomplished  by  successively  solving  an  equation 
for  one  of  the  unknowns  (not  P2)  and  substituting  the  resulting  expression  into  the 
remaining  equations.  The  system  is  thus  repeatedly  reduced  by  one  equation  and  one 
unknown  until  there  is  only  one  remaining  equation  with  P£  as  the  unknown. 

First  Equation  (C-l)  is  solved  for  V2. 


v2  -  tanB  (u2  -  U]P|/P2)  (C-5) 

The  following  Equations  are  obtained  by  substituting  the  expression  above  into 
Equations  (C-2)  through  (C-4) . 

P1  +  Piu3_2sin28  “  P2  +  8in28  u^p^/p  (C-6) 


2 

u^cos  B 


u2  -  sin  B  UjPj^Pj 


(C-7) 
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♦ 

4 


Ul2  "  U2  +  -  Zf)  "  tan2B  (U2  '  Vl/p2)2 


(C-8) 


Equation  (C-7)  is  solved  for  u2< 


U2  "  “l 


(cos2e  +  sin2e  P1/p2^ 


(C-9) 


The  system  is  reduced  to  two  equations  and  two  unknowns  by  substituting  Equation 
(C-9)  into  Equation  (C-8) .  Equation  (C-6)  is  unchanged  since  it  does  not  contain 


V 


2  2 
u^  sin  6 


-  1 


"1  \P1  p2/ 


Jl. 

Y 


(C-10) 


The  system  now  includes  Equation  (C-6)  and  (C-10) .  Equation  (C-6)  is  solved  for 
P2* 

p2  -  pi  +  PlUl2sin28  (1  - 
Substituting  Equation  (C-ll)  into  (C-10)  yields 

P22(piui2sin28  -  YP.jU^2sin2B  -  2yp1)  + 

p2(2yp1p1  +  2yP12u12sin2S)  -  P^u^sin^  -  YP^u^sin^  -  0. 

The  two  solutions  to  this  quadratic  equation  are 


(C-ll) 


(C-12) 


and 


P2  "  Pi 


P2  -  PX 


P1u12sin2B  (y+1) 

2yp1  +  P1u12sin2S  (y-1) 


(C-13) 


(C-14) 


Equation  (C-13)  Is  an  extraneous  solution.  Since  Mach  nunber  is  given  by 


m  2  2/  2  U1  P1 

M-  "  ui  /a-  “  • 


(C-15) 


Equation  (C-14)  may  be  written 
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Equation 


(Y+l)Mo2sin2g 

Pn  "  Pi  2  2  ’ 

(y-DM^  sin  6  +  2 

16)  is  the  desired  expression  for  p ,, . 


(C-16) 
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Appendix  D 

Derivation  of  the  Average  OIF  for  Random  Wavefront  Aberrations 
Consider  the  wave  phase  A(x,y)  given  by 

A(x,y)  -  Ac(x,y)  +  AR(x,y)  (D- 

where  A^  is  time-invariant  or  controlled  and  Ar  defines  a  spatially-stationary 
Gaussian  random  process.  The  purpose  of  this  appendix  is  to  find  an  expression 
for  the  expected  value  of  the  OTF  corresponding  to  the  wave  phase  A  given  by 
Equation  (D-l) . 

From  Equation  (2-28)  the  OTF  for  the  A  of  Equation  (D-l)  is  given  by 


^Af’Xf*"  I  //G0(x'*y')G0(x'+X>y'+y) 


x  exp{ik[Ac(x'+x,y'+y)  -  Ac(x',y')]}  (D-2 

x  exp{ik[AR(x'+x,y'+y)  -  AR(x',y')]}  dx'  dy'. 

Taking  the  expected  value  of  both  sides  of  Equation  (D-2)  yields 

oo 

E[T(yj,-^)]  m  Y  ff  GQ(x",y')G0(x'+x,y'+y)exp{ik[Ac(x'+x,y'+y)  -  Ac(x',y')]} 

mOO 

x  E^exp{ik[AR(x'+x,y'+y)  -  AR(x',y')  ])^dx'  dy'. 


Because  A^  is  spatially  stationary  the  expected  value  in  the  right  side  of  Equation 
(D-3)  is  a  function  of  neither  x'  nor  y'  and  may  be  extracted  from  the  integral. 

E[T(A?’Xf)^  "  Tc(xf *Xf)E(eXp{lk^R(X>X,y  +y)  "  AR(X  ,y  )^})  0>“4) 


CD 

TC(X?'Xf)  “  Z  If  C0(*',y')Go(x'+x,y'4y> 


x  exp{ik[Ac(x'+x,y'+y)  -  Ac(x',y')]}  dx'  dy' 

The  OTF  Is  the  OTF  which  would  be  obtained  if  A  were  completely  controlled.  To 
simplify  notation  the  following  terms  are  Introduced. 


U  -  AR(x',y')  V  -  AR(x'+x,y'+y) 


(D-6) 
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Equation  (D-4)  becomes 


‘  'c<w'ff)E(*lk<V’°>1-  <-»> 

U  and  V  are  dependent  random  variables  with  a  joint  probability  density  function 
given  by 


f(n,v) 


2ircr  a  /1-p  ‘ 
u  v 


exp  - 


where 


p  -  lE(UV)  -  P  u  ]/o  o 


U  V 


U  V 


(D-9) 


2  2 

and  where  y  ,  uv,  au  ,  and  av  are  respectively  the  expected  values  of  D  and  V  and 
the  variances  of  U  and  V.  The  bivariate  normal  density  function  given  by  Equation 
(D-8)  is  Implied  by  the  assumption  that  A^  defines  a  Guassian  random  process. 

The  numerator  of  the  right  side  of  Equation  (D-9)  is  the  convariance  function  <t> (x,y) 
Since  the  process  is  spatially-stationary,  *(x,y)  is  not  a  function  of  x'  or  y' , 
the  variances  of  U  and  V  are  equal,  and  the  expected  values  of  U  and  V  are  equal. 


2  2  2 
au  “  av  “  a 

2 

p  -  *(x,y)/a 

Without  loss  of  generality  it  is  assumed  the  mean  wave  phase  is  zero. 

y  ■  y  ■  0 

u  v 

Substituting  Equations  (D-10),  (D-ll),  and  (D-12)  into  Equation  (D-8)  yields 

, ,  1  T  cr^u^-lGuv+a^v^  1 

f  (u, v)  -  - -  exp - 7 — 5 - * 

L  2(o4-*2)  J 


(D-10) 

(D-ll) 

(D-12) 


(D-13) 


2ir  /a4-$‘ 

The  expected  value  in  the  right  side  of  Equation  (D-7)  may  now  be  expressed 


as 


E[eik(v-»)] .  n .  ;;eik(v-u)£(u  v)  du  dv 


(D-14) 


or 


EX 


2  2  2  2  4  2  -i 

a  u  -2$uTH-q  *v+2ik(g  -♦  )(v-u)|j^ 

2(o4-*2) 


(D-15) 
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4  2  2  2 

The  term  [$v  +  ik(o  -  4>) ]/ct  is  added  and  subtracted  from  the  numerator  to 
yield 


EX 


2i r/o4-$2  °° 


1 _  jfjfexpJ-  ^ ltr4-r d>v1-±lc (ct**— 1  (  dy 

.4  2-«  /  2  (a4-$2)  ’ 


4  ,2. 


4  ,2S,2,_2 


(D-16) 


x  exp 


_  q2v2+2ik(q4-$2)v-f  $vH.k(q4-$2) ]2/q2  j  dv 
2 (a4-*2) 


Notice  that  the  numerator  in  the  inner  integral  is  a  perfect  square.  With  the 
variable  substitution 


u  qu-fSv+ikfg4-1!12)  1/p 

w  -  /  L  2 

/2(o  -O 


(D-17) 


Equation  (P-16)  becomes 


EX  «  ==r-  /  e_W  dW 


(D-18) 


X  f  exp 


_  q^v^2ik(a4-t2)v-  [cnH-ik(o4-»2)  1 2/ q2  |  dv 
2(q4-*2) 


The  value  of  the  first  integral  is  /iT  and  the  exponential  term  in  the  second  integral 
may  be  simplified. 


EX 


1  ?  -[v2+2ik(q2-$)v]/2q2  -k2(q4-4>2)/2o2 

7Tva  J  e  dv  e 


(D-19) 


With  the  variable  substitution 


z  _  v+ik(q2-$) 


(D-20) 


Equation  (D-19)  becomes 
EX  - 


dZ  exp[- 


k2(q2-»)2-Pk2fq4-»2)1 
2c,2  ' 


(D-21) 


or 


EX  ■  erp[-k2(o2-*)] . 


(D-22) 
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Finally  Equations  (D-22)  and  (D-14)  are  substituted  into  Equation  (D-7). 

Equation  (D-23)  is  the  desired  expression  for  the  expected  value  of  the  OTF. 


(D-23) 
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Appendix  E 


Derivation  of  Equation  for  o 


The  covariance  function  and  variance  of  the  spatially-stationary  wave  phase  A 
are  given  respectively  by 

*(x,y)  ■  A(x',y')A(x'+x,y'+y)  (E-l) 


a2  -  *(0,0)  -  A2(x',y').  (E-2 

The  average  phase  is  defined  to  be  zero,  and  it  is  assumed  that  refractive  index 
changes  in  the  optical  path  are  small. 


A(x',y')  -  0 

n(x',y',z')  -  n  -  «  1 


The  above  conditions  imply  that  A  is  given  by 


Mx'.y')  -  / u(x',y',z')  dz'. 
0 


(E-5) 


In  this  appendix  expressions  are  derived  for  4>(x,y)  and  o  .  It  is  assumed  that 
the  refractive  index  variation  u  has  a  Gaussian  covariance  function. 


u(x',y',z')w(x'+x,y'+y,z'+z)  «  v  exp 


From  Equations  (E-l)  and  (E-5)  the  covariance  function  *  is  given  by 


L  L 

*(x,y)  -  J y(x',y',z')dz' J  u(x'+x,y'+y,z2)dz2 
0  0 


(E-7) 


$(x,y)  ■  fj  w(x',y',z')y(x'+x,y'+y,z2)dz2dz> 
00 


With  the  variable  substitution 


z  ■  z2  -  z 
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Equation  (E-8)  becomes 


L  L-z' 


i(x,y)  -J*  f  y  (x',y',  z')y(x'+x,y '+y,z'+z)  dz  dz'. 
0  -z 


(E-10) 


Substituting  Equation  (E-6)  into  Equation  (E-10)  yields 


$(x,y)  -  y  exp 


2  2 

exp(-z  / az  )  dz  dz' 


(E-ll) 


2  2 

4>(x,y)  -  la  y  exp 
z 


where  I  is  defined  by 


(E-12) 


L  L-z' 


I  ■  —^2  f  f  exp(-z2/az2)  dz  dz'. 
az  0  -z 


With  the  normalized  variables 


u  ■  z/a  v  ■  z'/a  w  -  L/a 

z  z  z 


Equation  (E-13)  becomes 


w  w-v  2 
I  mJ  J*  e  U  du  dv. 
0  -v 


(E-13) 


(E-14) 


(E-15) 


To  evaluate  I,  Equation  (E-15)  is  expressed  in  terms  of  the  error  function,  erf, 
and  the  error  function  complement,  erfc.  These  functions  are  given  by 


erf(v)  -  /e“u  du 

0  -v 


erfc(v)  -  1  -  erf(v)  -  -jL  Je^du  »  j  e-u  c 


erf(O)  ■  1  -  erfc(»)  ■  0 


(E-16) 


(E-17) 


(E-18) 
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erf(®)  ■  erfc(O)  -  1 

The  inner  integral  in  Equation  (E-15)  may  be  written 

2  “  .2  ~Y  ..2 

|  e  U  du  -  e  * 

-v  -®  0  w-v 


J  e~U  du  *  Je  U  du  +  Je  U  du  -  J  e  U  du  -  Je  U  du 


■  ~  [2  -  erfc  (w-v)  -  erf(v)] 
Substituting  Equation  (E-20)  into  (E-15)  yields 


I  -  2w- J  erfc(w-v)  dv  -  J*  erfc(v)  dv 

0  0 

•  A“  [w- y* erfc(v)  dvj. 

The  error  function  complement  has  the  property^® 


w 


(E-19) 


(E-20) 


(E-21) 


00  «  2 
^erfc(v)dv  ■  ■^^’(v-w)e 


dv. 


(E-22) 


Equation  (E-21)  becomes 


I 


Aw 


-  2/ve"v  dv 
0 


■  Aw 


1  +  e 


-w 


dv 


(E-23) 


2 

-  Aw  erf(w)  -  1  +  e-w  . 

Substituting  Equations  (E-14)  and  (E-23)  into  Equation  (E-12)  yields  the  general 
expression  for  4. 


*(x,y) 


-  [a'  (L/a8)erf(L/a8)  -  1  +  e”(L/*z)  ]az2y2exp  ~(^~2  +  JL2^ 


(E-24) 
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If  a  «  L,  then  Equation  (E-24)  becomes 
z 


$(x,y)  •  j*  h  u  exp 


[tr$l 


*  y 

Finally,  If  u  is  isotropic,  then  the  a's  are  all  equal, 

a  ■  a  ■  a  “a, 
x  y  z  ’ 

and  the  covariance  function  and  variance  of  4  are  given  respectively  by 

$(x,y)  -  /T  a  L  u2  exp[-(x2  +  y2)/a2] 


and 


2  r~  2 

a  ■  n  a  L  p  . 


(E-25) 


(£-26) 


(E-27) 


(E-28) 
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Appendix  F 

Derivation  of  £  f°r  the  S tine-Winovich  Experiment 

The  purpose  of  this  appendix  is  to  find  the  average  theoretical  ratio  Q 
of  the  radiant  power  emerging  from  the  collection  aperture  to  that  emerging  from 
the  source  aperture  for  the  Stine-Winovich  experiment. 

From  the  problem  geometry  and  the  assumption  of  an  isostropic  boundary  layer 
disturbance  both  the  radiant  emlttance  of  the  source  aperture  0(r)  and  the  image 
irradiance  I(R)  are  rotationally  synmetric.  The  ratio  ^  is  therefore  given  by 


Q 


2ir  /  I(R)  R  dR 
0 _ 

*a 

2*  f  0(r)  r  dr 
0 


(F-l) 


It  is  assumed  that  0(r)  is  constant  over  the  source  aperture,  and  for  simplicity 
the  total  radiant  power  through  the  source  aperture  is  taken  to  be  unity.  0(r) 
and  Q  are,  therefore,  given  respectively  by 


and 


0(r) 


0  <  r  <  ra 


0  otherwise 


Q 


1 (R)  R  dR. 


In  the  spatial-frequency  domain  I(R)  and  0(r)  are  related  by 


(F-2) 


(F-3) 


S,  (v)  -  S  (v) t (v)  (F-4) 

1  o 

where  and  Sc  are  the  Fourier  transforms  of  I  and  0  respectively  and  t  is  the 
system  OTF.  Since  the  system  is  rotationally  symmetric,  S. ,  S  ,  and  t  are  all  real. 
The  Fourier  transform  pairs  I,  ,  and  0,  S0  have  the  followln§  relationships. 

00 

1(1)  -  21 rf  S1(v)J0(2irvR)  v  dv  (F-5) 


a 
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Sq(v)  ■  2ir  f  0(r)JQ(2irvr)  r  dr 


(F-6) 


f  V2™0  r 


dr 


The  following  equation  is  obtained  by  substituting  Equations  (F-4) ,  (F-5) ,  and 
(F-6)  Into  Equation  (F-3) . 


2  rb  oo  ra  __ 

Q  -  f  f  J  JQ(2irvr)  r  dr  t(v)Jq(2itvR)  v  dv  R  dR 


r  0  0  0 
a 


(F-7) 


By  changing  the  order  of  Integration,  Equation  (F-7)  may  be  written 

.  2  •  rb  Ja 

Jg(2ifvR)  R  dR  J  Jg(2irvr)  r  dr  r(v)  v  dv 

“  0 


ra  0  0 


Bessel  functions  have  the  property 


27 


x  dl  [xJl(x)1  “  J0(x)- 


(F-8) 


(F-9) 


In  Integral  form  Equation  (F-9)  may  be  written 

x 


f  rJQ(r)dr  -  xJ^x), 


(F-10) 


With  the  formula  above.  Equation  (F-8)  becomes 

2r  00  — 

Q  -  7 —f  ^  J1(2irvrb)  J  (2Trvr  )  dv. 
b  0 

The  upper  limit  of  the  integral  may  be  changed  to  D/lf  since e 

t(v)  ■  T(ff)  ■  0,  r>D 


(F-ll) 


(F-12) 


where  D  Is  the  telescope  diameter. 
With  the  variable  substitution 


u  ■  2irvr. 


(F-13) 
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Equation  (F-ll)  becomes 

2nDr. 

- b  ra 

« ■  /  **• 

u  _a 

rb 


(F-14) 


The  average_OTF  t  is  now  examined.  In  terms  of  its  random  and  diffraction-limit 
components  t  is  given  by 


T(Xf)  “  TR2(Xf^T0(Xf)  (F_15> 

where  a  "2"  is  included  in  the  subscript  of  tr  to  indicate  that  two  TBL's  are 
crossed  by  the  optical  beams.  It  is  assumed  that  the  two  TBL's  are  independent 
but  statistically  identical.  The  random  optical  phase  for  a  beam  crossing  both 
TBL's  is  given  by 


-  V  +  V  <r-16) 

where  Aj^  and  ARb  are  the  random  phases  caused  by  the  first  and  second  TBL's  respec¬ 
tively.  Since  Ar^  and  ARb  are  independent  Gaussian  random  variables,  the  variance 
of  Ar2  is  equal  to  the  sum  of  the  variances  of  Aj^  and  ARb.  If  the  variances  of 
ARa  and  ARb  are  equal  to  o  ,  then  the  variance  of  AR2  is  given  by 

2  2.2  „  2 

o2  -  o  +  a  -  2o  (F-17) 


The  equation  for  tr2  is  therefore 


(F-18) 


where  tr  is  the  random  OTF  for  a  beam  crossing  only  one  TBL.  Substituting  Equation 
(F-15)  and  (F-18)  into  Equation  (F-14)  leads  to  the  final  expression  for  T?. 


2wDrt 

~Tf~ 

Q  “  / 

o 


2  2J1<“  t' 


du 


(F-19) 
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Appendix  G 

Derivation  of  the  Radiant  Output  Power  of  a  Shearing  Interferometer 

as  £  Function  of  the  MTF 

Figure  58  shows  the  superimposed  beams  for  a  shearing  Interferometer.  It  Is 
assumed  that  the  beams  are  monochromatic  and  perfectly  coherent.  The  beam  I  is 
displaced  from  the  beam  Z^  by  (x,y) ,  and  the  relative  phase  of  the  two  beams  is  A 
In  general  A  is  a  function  of  either  (x,y)  or  time  t  or  both  depending  on  the  inter¬ 
ferometer  design. 

The  pupil  functions  for  Z  and  Zj  are  given  respectively  by 


*',y'>  -{ 


AeikA(x'.y'>, 
0  otherwise 


_  .  .  ..  +  «•  +  »>  +  M.  (x',y-)c£l  (M 

V  ,y  *  (0  otherwise. 

Since  the  condition  (x'.y'JeZ^  implies  the  condition  (x'  +  x,y'  +  y)eZ,  the  pupil 
function  Gi  may  be  written 


G.Cx'.y')  -  G(x'  +  x,  y'  +  y)e 


(G-3) 


The  pupil  function  for  the  interferometer  output  Gj  is  the  sum  of  G  and  G^, 

G2(x',y')  -  G(x'.y')  4-  G^x'.y') 


G2(x%y')  -  G(x',y')  +  G(x'  +  x,y'  +  y)e““. 

The  radiant  power  density  of  the  interferometer  output  is 

I(x',y')  -  G2*<x',y')G2(x',y') 

From  Equation  (G-5)  Equation  (G-6)  becomes 

I(x',y")  -  G*(x',y')G(x',y')  +  G*(x'+x,y'+y)G(x'+x,y'+y> 

+  G*(x',y')G<x'+x,y>y)eikA+  G*(x'+x,y'+y)G(x',y')e“ikA. 
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V' 


FIGURE  58 

CROSS  SECTION  OF  SUPERIMPOSED  BEAMS  FOR  A  SHEARING  INTERFEROMETER 
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pew  -  2EA*A|l  +  |T<Y£,-^)|cos[n(x,y)+kA]}.  (G-12) 

In  a  real  system  finite  bandwidth.  Inequality  of  power  of  the  two  beams,  and 
partial  Incoherence  of  the  beams  prevent  Equation  (G-12)  from  being  realized.  It 
can  be  shown,  however,  that  even  with  the  error  sources  mentioned  above,  pow  Is 
still  directly  proportional  to  the  amplitude  of  the  cosine  term.  Equation  (G-12) 

Is  modified  to 


pow  -  a  +  8 1  T  (xp'xf  ) I  cos  £ (*,y)  +  kAl  (G-13) 

where  a  and  8  are  constants  (8<a) ,  and  is  the  phase  change  caused  v /  partial 
incoherence. 

For  the  comer-cube  design  the  phase  term  kA  is  directly  proportional  to  time. 
With  the  x-axis  as  the  direction  of  shear.  Equation  (G-13)  becomes 

pow  -  a  +  6|t(yL,0)|cos[$(x)  +  w  t]  (G-14) 

where 

*(x)  -  ♦1  +  n(x,0)  (G-15) 

and  ui  is  the  angular  frequency  of  the  relative  phase  of  the  two  beams. 
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